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1.0  INTRODUCTION 

Much  of  the  testing  and  research  at  Arnold  Engineering  Development 
Center  (AEDC)  evaluates  quantitatively  the  flows  characteristic  of  re¬ 
active  propulsion  systems.  Although  the  range  of  conditions  experienced 
in  such  studies  is  diverse,  extremely  high  velocities,  turbulence  levels, 
and  temperatures  are  often  encountered  in  a  single,  chemically  reacting 
medium.  In  recent  years  attention  has  been  given,  therefore,  to  nonin- 
trusive  instrumentation  systems  which  can  determine  the  velocity, 
density,  temperature,  or  species  concentration  in  such  flows  without 
contact  with  the  fluid  medium. 

Particular  attention  has  been  given  to  techniques  for  determining 
the  nature  of  the  velocity  field  in  turbulent  shear  flows  since  con¬ 
vection  mechanisms  contribute  to  the  development  of  the  temperature 
field  and  chemical  reactions  and  to  the  thrust  of  a  propulsive  device. 

The  predominant  nonintrusive  instrument  for  this  purpose,  the  laser 
velocimeter  (LV) ,  has  been  used  to  measure  (1)  the  mean  velocity  and 
turbulence  structure  in  heated  and  unheated  jets  at  velocities  rang¬ 
ing  from  the  low-subsonic  to  near-hypersonic  regimes  (Refs.  I  and  2), 

(2)  the  recirculation  and  mixing  in  reacting,  coaxial  hydrogen-air  jets 
(Refs.  3  and  4),  (3)  the  fully  developed  turbulent  pipe  flows  (Ref.  1), 
and  (4)  the  change  in  centerline  exhaust  velocity  of  a  large  turbofan 
engine  with  throttle  setting  (Ref.  51. 

Although  each  of  these  experiments  was  successful  in  meeting 
its  stated  objectives,  a  critical  assessment  of  various  factors 
influencing  the  results  revealed  many  potential  sources  of  error 
in  the  measurements.  Specifically,  the  accuracy  of  the  data  may  be 
influenced  by  (1)  the  LV  design,  (2)  experimental  procedures,  (3) 
data  acquisition  and  reduction  techniques,  (4)  Instrumentation  un¬ 
certainties,  and  (5)  particle  dynamics  effects.  Detailed  discussions 
of  these  possible  errors  are  given  in  several  AEDC  publications  (Refs. 
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1,  2,  and  6  through  10)  and  in  the  open  literature  (Refs.  11  through 
14).  The  present  study,  however,  concerns  the  elimination  of  two 
sources  of  possible  serious  errors — statistical  bias  and  large  par¬ 
ticle  lags.  Since  the  overall  LV  development  at  AEDC  includes  de¬ 
velopment  of  burst  processors  (Refs.  15  and  16)  as  well  as  LV  utili¬ 
zation,  the  topics  of  statistical  bias  and  large-particle  discimination 
techniques  apply  to  the  entire  scope  of  laser  velocimetry  at  AEDC, 

1.1  BACKGROUND 

A  dual-scatter  laser  velocimeter,  such  as  shown  in  Fig.  1,  consists 

of  (1)  a  laser,  which  provides  a  monochromatic,  coherent  light  source, 

(2)  illuminating  optics,  which  split  the  laser  beam  into  two  or  more 

essentially  equal  intensity  beams  then  refocus  them  to  form  an  optical 

probe  volume,  (3)  collector  optics,  which  collect  light  scattered  by 

particles  traversing  the  interference  fringes  established  in  the  probe 

volume,  and  (4)  system  electronics,  which  amplify  the  signals  and 

determine  the  frequency,  f^,  at  which  the  particles  intercept  the  fringes. 

Since  the  fringe  spacing,  K. ,  is  determined  solely  by  the  illuminating 

5 

optics  (see  Ref.  17),  the  particle  velocity  normal  to  the  plane  of  the 
fringes,  V,  is  simply  the  product  K^f^.  If  the  particles  are  small 
enough  to  nearly  duplicate  the  spatial  and  temporal  variations  of  the 
fluid  velocity,  an  extremely  accurate  measurement  of  local  fluid  velocity 
can  be  obtained. 

When  the  concentration  of  particles  large  enough  to  be  detected  by 
the  LV  is  low,  as  in  most  gaseous  flows,  data  will  be  obtained  in¬ 
termittently  as  single  particles  are  swept  through  the  probe  volume.  In 
studying  the  sampling  characteristics  of  laser  velocimete'rs  operated 
in  such  an  individual  realization  mode,  Mayo  (Ref.  18)  observed  that  the 
interval  between  particle  arrivals  was  described  by  a  Poisson  prob¬ 
ability  distribution  having  a  rate  parameter  which  was  proportional 
to  the  flux  of  particles  passing  through  the  probe  volume.  He  con- 
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Figure  1.  Dual-scatter  laser  velocimeter. 
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eluded  that,  since  the  sampling  rate  and  measured  velocity  were  inter¬ 
related,  individual  realization  laser  velocimeters  (IRLV's)  would 
overestimate  the  mean  velocity  at  a  point.  A  similar  conclusion  was 
independently  reached  by  McLaughlin  and  Tiederman  (Ref.  19),  who 
termed  the  phenomenon  statistical  bias.  Barnett  and  Bentley  (Ref.  20) 
showed  that  the  mean  velocity  was  biased  in  direct  proportion  to  the 
turbulence  intensity,  thus  identifying  a  potentially  serious  source 
of  error  in  shear  flows. 

Since  the  publication  of  these  pioneering  studies,  considerable 
research  has  been  devoted  to  determining  the  conditions  under  which 
bias  occurs  and  to  developing  data  acquisition  or  processing  techniques 
to  eliminate  bias  (Refs.  21  through  25).  Considerably  less  work  has 
been  devoted  to  the  experimental  verification  of  the  magnitude  of  bias 
(Refs.  26  and  27).  Data  obtained  at  AEDC  (Ref.  1)  have,  moreover,  not 
shown  any  conclusive  evidence  that  the  theoretically  predicted  bias 
exists.  Despite  this  lack  of  experimental  verification  many 
experimentalists  routinely  correct  IRLV  data  for  bias  and  some 
commercial  burst  counters  developed  for  LV  applications  (Ref,  28) 
process  IRLV  data  using  algorithms  developed  to  remove  bias. 

The  magnitude  and  existence  of  statistical  bias  are  of  par¬ 
ticular  concern  to  the  LV  development  activities  at  AEDC  since  (1) 
data  obtained  in  turbulent  and  transient  flows  could  be  especially 
influenced  by  bias  and  (2)  burst  processor  development  activities 
should  consider  direct  removal  of  any  significant  data  biases. 

In  using  an  IRLV,  a  data  set  containing  measurements  of  the 
velocity  of  many  individual  particles  is  obtained.  The  response 
of  each  particle  to  the  spatial  or  temporal  velocity  variations 
in  the  fluid  is  dictated  by  the  trajectory,  shape,  size,  and  density 
of  the  particle.  If  these  parameters  could  be  controlled  for  all 
particles  in  a  data  set,  it  would  be  possible  to  accurately  predict 
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how  closely  the  LV  measurements  correspond  to  the  actual  fluid 
velocity  at  a  point.  In  applications  where  the  flow  is  seeded  with 
a  single  size  and  type  of  particle,  this  condition  may  be  approached 
but  not  reached  because  of  particle  agglomeration  or  the  presence 
of  other  particles  in  the  fluid.  Other  applications  preclude  active 
seeding  of  the  flow  so  that  natural  particulates  of  widely  varying 
properties  contribute  to  a  data  set. 

Ideally,  then,  an  IRLV  should  measure  particle  size  as  well  as 
velocity  if  the  best  accuracy  of  the  instrument  is  realized.  Sig¬ 
nificant  advances  in  particle  size  determination  by  laser  velocimetry 
have  been  made  by  Farmer  (Ref.  29)  and  later  by  Roberds  (Ref.  30).  In 
both  studies,  the  signal  visibility  of  the  individual  Doppler  burst  is 
determined  and  related  to  the  particle  size  through  scattering  theory. 
Under  carefully  controlled  conditions  (Ref.  29),  the  visibility  is  in¬ 
dependent  of  index  of  refraction  and  the  dimension  of  the  particle 
normal  to  the  probe  volume  fringe  planes  can  be  determined.  Such  results, 
however,  are  not  unique  if  the  visibility  is  low,  i.e.,  a  single 
visibility  value  is  characteristic  of  several  particle  sizes.  The 
visibility  method  is  at  present  additionally  limited  to  relatively 
low  frequencies  (10  MHz  or  less),  making  it  of  limited  utility  since 
high  velocities  or  acousto-optical  modulation  (Ref.  31)  of  the  signal 
to  determine  flow  direction  can  lead  to  much  higher  signal  frequencies. 
Recently,  moreover,  several  Investigations  (Refs.  32,  33,  and  34)  have 
raised  serious  questions  as  to  the  validity  of  the  visibility  concept 
in  backscatter  applications. 


Because  of  these  limitations,  use  of  signal  visibility  to  deter¬ 
mine  particle  size  is  not  presently  suitable  for  many  applications. 
If,  however,  the  criterion  of  simultaneous  particle  size  and  velocity 
determination  is  relaxed,  an  alternate  approach  may  be  considered. 

Reference  2  presents  LV  data  obtained  from  a  highly  accelerated 
A  2 

(1.3  x  10  ft/sec  ),  naturally  seeded  flow.  In  many  cases  there  were 
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such  wide  variations  in  measured  velocities  in  each  data  set  that  the 
mean  velocity  was  useless  as  an  indicator  of  local  fluid  velocity. 

Note,  however,  that  the  upper  velocity  limit  of  the  measured  probability 
distributions  approached  the  gas  velocity  in  acceleration  regions, 
whereas  the  lower  limit  was  near  the  gas  velocity  after  a  sudden  decel¬ 
eration.  Since  these  results  are  typical  of  the  particle  dynamic 
behavior  of  an  ensemble  of  particles  having  a  wide  distribution  of  sizes 
and  shapes,  an  appreciable  increase  in  LV  accuracy  may  result  if  signals 
arising  from  large  particles  can  be  excluded  from  a  data  set. 

Andrews  and  Seifert  (Ref.  35)  attempted  to  size  particles  by  the 
amplitude  of  light  scattered  by  a  reference-beam  LV.  Pfeiffer  and 
von  Stein  (Ref.  36)  reported  that  discriminating  against  signals  whose 
amplitude  exceeded  a  predetermined  threshold  value  led  to  significant 
improvements  in  dual-scatter  LV  accuracy  in  accelerated  flows.  Again, 
the  instrument  manufacturers  (Ref.  28)  seized  on  these  results  and  have 
produced  data  processors  that  incorporate  signal  amplitude  discriminators. 
Unfortunately,  these  studies  do  not  provide  the  detail  necessary  for,  ' 
a  p>u.0fu.,  accurately  estimating  the  particle  sizes  being  excluded  from 
consideration  in  obtaining  data.  If  the  threshold  is  set  too  low,  data 
acquisition  rates  are  adversely  affected,  and  too  high  a  setting  will 
allow  too  many  large  particles  to  be  considered.  Moreover,  the  com¬ 
plexity  of  the  Mie  scattering  process  is  such  that  a  slight  variation 
in' either  the  particle  characteristics  or  the  LV  system  could  sig¬ 
nificantly  alter  the  results.  Although  the  concept  of  discrimination 
on  the  basis  of  signal  amplitude  is  intuitively  and  operationally 
attractive,  further  study  is  warranted  before  incorporating  this  con¬ 
cept  into  future  LV  systems  and  processors. 

1.2  SCOPE  OF  THE  INVESTIGATION 

The  present  analytical  and  experimental  study  was  undertaken  to 
determine  (1)  the  magnitude  of  statistical  bias  in  individual  realization 
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laser,  velocimeter  data  and  (2)  the  suitability  of  discriminating  against 
signals  from  large  particles  on  the  basis  of  signal  amplitude.  Because 
of  the  importance  of  these  concepts  in  obtaining  improved  measurements 
in  the  shear  flows  prevalent  in  many  tests,  the  experimental  data 
were  obtained  in  the  shear  layer  formed  by  a  subsonic  jet  exhausting 
into  a  quiescent  medium.  To  minimize  the  uncertainty  of  particle  dy¬ 
namic  effects  on  the  data,  all  particulates  exceeding  0.1  pm  in  diameter 
were  filtered  from  the  air  supply  and  the  flow  then  seeded  with  particles 
having  known  characteristics  by  either  a  fluidized  bed  seeder  or  aerosol 
generator.  Entrainment  of  particles  from  the  ambient  air  was  reduced 
by  confining  the  developing  jet  flow  field  within  an  enclosure  having 
transparent  walls  to  permit  LV  utilization.  The  length  of  the  enclo¬ 
sure  was  chosen  to  enhance  recirculation  along  the  edge  of  the  jet, 
thus  increasing  uniformity  of  seed  concentration  across  the  jet. 

In  Section  2.0,  pitot-static,  hot-wire,  and  one-component  LV  data 
of  the  axial  and  radial  variation  of  velocity  and  turbulence  intensity 
in  the  jet  near  field  are  presented  for  axial  location  to  nozzle  exit 
diameter  ratios  of  2,  A,  6,  and  8.  The  data  are  compared  to  statistical 
bias  theories  and  a  bias  simulation  program  developed  for  the  studies 
to  determine  (1)  the  magnitude  of  statistical  bias,  (2)  its  relationship 
to  seeding  rate  and  turbulence  structure,  and  (3)  the  factors  influential 
in  the  reduction  of  bias  in  jet  flows. 

Section  3.0  summarizes  those  conditions  that  influence  the  amplitude 
of  signals  detected  by  dual-scatter  laser  velocimeters .  Specifically, 
the  influence  of  particle  size,  index  of  refraction,  illuminating 
intensity  variation,  and  LV  collector  optics  on  the  amplitude  of  the 
Mie  scattering  intensity  at  the  phototube  is  discussed.  Photomultiplier 
characteristics  and  other  signal-amplification  electronic  components 
are  then  considered  in  arriving  at  the  signal  voltage  level  observed 
by  the  data  processor.  The  theoretical  signals  are  compared  to  those 
observed  from  the  subsonic  jet  for  particles  of  varying  sizes.  An 
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assessment  of  the  potential  for  discriminating  against  large  particles  - 
by  signal  amplitude  discrimination  is  made. 

2.0  STATISTICAL  BIAS  STUDY 


To  experimentally  evaluate  the  magnitude  of  the  statistical  bias 
presumed  present  in  IRLV  measurements  of  the  mean  velocity  and  higher 
velocity  moments,  those  flow  conditions  that  lead  to  biased  results 
must  be  identified  and  controlled.  The  bias  observed  in  a  given  application 
depends  upon  the  turbulence  level,  data  acquisition  system,  data  re¬ 
duction  procedure,  and  concentration  of  the  particles  from  which  the 
data  are  obtained.  Moreover,  bias  can  only  be  determined  relative  to 
other  instrumentation  systems.  Therefore,  the  absolute  accuracy  of 
the  LV  and  supporting  instrumentation  must  be  within  certain  bounds  if 
the  effect  is  to  be  observed. 


2.1  THEORETICAL  CONSIDERATIONS 


The  bias  model  proposed  by  Barnett  and  Bentley  (Ref.  20)  has  been 
extended  to  identify  the  primary  bias-creating  flow  conditions  and  the 
anticipated  magnitude  of  the  bias.  This  model  was  chosen  because  it 
leads  directly  to  a  quantitative  assessment  of  the  statistical  bias  and 
it  directly  relates  arithmetic  averages  to  the  "unbiased11  temporal 
averages  which  form  the  bias  of  many  turbulent  theories. 


The  temporal  average  of  any  moment,  f,  of  the  velocity  may  be 
obtained  over  an  interval  T  from 


f 


I 

T 


t  -*-T 
0 

/  f{t)  dt 

t 

o 


(1) 
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IRLV  systems,  however,  measure  the  velocity  only  at  those  discrete  times 

A 

when  a  particle  is  crossing  the  LV  probe  volume,  so  that  f  can  only  be 
approximated  by  a  numerical  integration  such  as 


f  = 


M 

L  r,At, 

1=  j 
M 

1-  i 


(2) 


In  Eq.  (2)  the  f  are  velocity  quantities  evaluated  for  each  of  M  discrete 
LV  measurements  and  the  At^  are  the  times  separating  the  respective 
measurements.  When  the  intervals  between  successive  measurements  are 
small,  f  will  be  a  good  approximation  of  f.  The  numerical  integration 
will*  furthermore,  reduce  to  the  arithmetic  average  of  the  IRLV  data 
set  under  certain  conditions.  For  example,  if  the  interval  between 
successive  samples  is  constant,  Eq,  (2)  will  reduce  to 


f 


M 

=  1  T 

M  " 
1=  1 


fi  5  <  f 


(3) 


where  the  angle  brackets  denote  the  arithmetic  average  of 


f 


i* 


The  nature  of  statistical  bias  in  individual  realization  ve- 
locimetry  is  also  evident  from  Eq.  (2)  since  a  Reynolds  decomposition 
of  the  velocity  moment  and  time  interval  into 


f j  =  <  f  >  + 


At  j  -  <  At  >  +  At; 


allows  the  numerical  approximation  of  the  temporal  average  to  b< 
written 

r  -  <t>  [i  +  h[4i] 

where  the  correlation  function  is  defined  as 

m  , ,  ,  , 

RfA  =  I  V  f.  At 

M  <  f  >  <  At> 


(4) 

(5) 

(6) 

(7) 
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It  is  apparent,  then,  from  Eq.  (6)  that  if  a  correlation  exists  between 

the  velocity  and  the  interarrival  time  between  particles,  the  arithmetic 

average  will  be  biased  with  respect  to  the  approximate  temporal  average. 

On  the  other  hand,  if  the  velocity  moment  and  time  between  samples,  f ^ 

and  At  ,  are  statistically  independent,  as  in  the  case  where  At  is 

i 

constant,  the  two  averages  are  identical. 

The  relationship  of  the  velocity,  to  the  time  interval,  At^  in 

a  time-varying  flow  can  be  established  under  certain  conditions.  Ac¬ 
cordingly,  assume  that  (1)  the  particle  number  density  (or  concentration), 
N,  is  constant  throughout  the  flow  field,  (2)  the  LV  probe  volume  is 
spherical  so  that  it  presents  a  constant  cross-sectional  area.  A,  to  all 
particles  traversing  it,  (3)  the  fluid  density,  p,  is  constant,  and  (4‘) 
all  particles  are  travelling  at  the  fluid  velocity  (negligible  lag). 

Then  the  number  of  particles  that  enters  the  probe  volume  in  At^  is 

i  +At 

SN-  -  NA  J  V (i)  dt  =  MAt- V. 

i  '  ‘  (8) 

O 

where  is  the  mean  fluid  velocity  during  the  time  between  successive 

samples,  At^  It  is  implicitly  assumed  in  Eq.  (8)  that  fluid  is  not 
recirculated  through  the  probe  volume  so  that  each  particle  can  enter 
the  probe  volume  only  once.  If  every  particle  velocity  is  processed 
by  the  LV  (i.e.,  6N±  —  1 ) s  then  the  interarrival  time  and  the  mean  fluid 
velocity  during  At^  are  related  by 

Att  =  _ 1 _ 

N A  V  (9) 

1 

If  At  is  small,  then  V(t)  is  essentially  constant  over  the 
interval  and  V_^  =  the  measured  particle  velocity,  and  from  the 
discussion  following  Eq.  (7)  it  is  expected  that  the  LV  data  set  would 
be  biased.  For  a  time-varying  flow,  characterized  by  a  single  frequency, 
this  condition  is  readily  obtained  if  At^  is  much  less  than  the  period 
of  the  flow  oscillation.  Turbulent  flows,  however,  are  typified  by  an 
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energy-cascading  process  extending  over  a  spectrum  of  frequencies. 
Accordingly,  it  is  more  difficult  to  determine  whether  a  relationship 
exists  between  data  rate  and  the  temporal  velocity  variation  in  turbulent 
flows,  although  such  a  relation  clearly  exists  for  simple  flow  pulsations. 
Some  insight  is  gained  if  the  variance  of  the  flow  velocities  measured  over 
At^  is  compared  to  the  variance  of  the  instantaneous  velocities.  In 
Appendix  A  these  variances  are  related  by 


tv, 


vy 


Ai 


f  R  (r)  r 

1  -  1  dr 

J  L 

At  J 

(10) 


where  the  autocorrelation  coefficient  is  defined  as 


aw  -  [(v,  -  v)2]  1  f  (V[t)  -  v)(vu  +  r)  -  v) 


dt 


(11) 


When  all  the  At^  time  intervals  are  sufficiently  small  R(t)  approaches  unity 
for  all  t  much  less  than  At  (Ref.  37).  Then  it  follows  immediately  from 


A  AH  A  rt  ' 

Eq.  (10)  that  (V^-V)  -  (V^-V)  .  An  obvious  condition  leading  to  this 

equality  is  for  the  mean  velocity  over  Ati  to  equal  the  measured  velocity 

(Vi  a  V.)  as  noted  earlier.  Conversely,  if  all  At^  are  large  enough  that 
R(t)  approaches  zero  faster  than  t  approaches  At^,  Eq.  (10)  reduces  to 

A^n  2 1 

(V.  -  V)2  (V,  -  V)2 

At.  (12) 


where  the  integral  time  scale  is  defined  by 


Lt  = 


lim 


T  =  Th 


/  R 


(r)  dr 


03) 


The  time  intervals  are  much  larger  than  the  integral  time  scale  then,  so 


A  A 


from  Eq.  (12),  the  quantity  (V^-V)  approaches  zero  and,  from  the  theorem 
of  the  mean  (Ref.  37),  the  mean  velocities  averaged  over  the  interarrival 
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times  approach  the  temporal  mean  velocity.  For  this  case  Eq.  (9)  shows 
that  the  particle  interarrival  times  are  constant  so  that  no  bias  of  the 
velocity  moments  computed  from  IRLV  data  would  be  expected.  The  data, 
rates  are,  however,  too  low  for  practical  use  of  a  laser  velocimeter,  as 
is  shown  below. 


The  magnitude  of  the  expected  bias  may  be  estimated  by  considering 
the  short  interarrival  time  case  which  requires  that  Eq.  (9)  become 

Ai  =  1  — 

'  N  A  V 


(14) 


and,  from  Eq.  (2), 


M 

Evv, 

f  =  j=1 _ 


M 

El/V 

i- 1 


(15) 


Note  that  for  small  Ap  and  a  sufficiently  large  data  set  f  =  f  is  an 
accurate  approximation  so  that  the  "true"  temporal  average  and  the 
finite-difference  integration  are  equivalent.  In  the  following,  such 
equivalence  is  assumed.  Therefore,  the  relation  for  p  =  is 


f 


V  =  _J_ 

<±> 

v 


(16) 


so  that  the  fluid  velocity  is  obtained  from  the  harmonic  average  of  the 
individual  velocity  measurements: 

V  =  =  _ ! _ 

<I>  M 

V 

The  more  general  case  of  the  bias  of  the  r1"11  central  moment  of 
velocity  is  readily  obtained  with  Eq.  (14)  for  short  interarrival  times 
if  it  is  noted  from  the  binomial  theorem  that 

r,  =  (v,  -  v)r  £(-i)s  (0  vp  Vs  <18) 

S”  o 


Et^ 


1=1 


v . 


(17) 
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where  (]£)  denotes  the  binomial  coefficients.  From  Eqs.  (15),  (17),  and 
(18),  it  follows  that 


(V  -  VV  =  Y,  (-n'h  (*>  <  v^1  >  V 


-+1 


(19) 


For  the  second  central  moment  r  =  2  and  Eq.  (19)  shows  that  arithmetic 
averages  of  the  mean  velocity  are  biased  with  respect  to  the  temporal 
averages  since 


<  V>  =  V  (1  +  4) 


(20) 


where  the  turbulence  intensity  is  defined  as 


Oxr 


(21) 


Similarly,  evaluating  Eq.  (19)  for  the  third  central  velocity  moment 
(r  =  3)  shows 


'<v> 


14  °\ 


(1  -  oi 


+  Sff  ~)  2 


(22) 


where  the  skew  coefficient  is  defined  by 


(23) 


Consequently,  the  bias  of  any  moment  computed  by  arithmetic  averages 
is  proportional  to  the  next  higher  moment  of  the  velocity. 


The  estimated  statistical  bias  of  the  mean  velocity,  <V>,  and 
turbulence  intensity,  °<v>»  when  the  average  interval  between 
realizations  is  much  less  than  the  flow  integral  time  scale,  is  shown 
in  Fig.  2.  The  arithmetic  average  velocity  always  overpredicts  the 
temporal  mean  velocity,  but  the  arithmetic  turbulent  intensity  is 
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strongly  affected  by  the  skewness  of  the  velocity  probability  distribution. 
Note  that  for  a  Gaussian  distribution  of  velocities  (S  =0),  local 
turbulent  intensities  computed  from  arithmetic  averages  will  under¬ 
estimate  the  true  value. 


The  relationships  derived  above  are  valid  only  for  a  one-dimensional 
flow  or  an  LV  whose  fringe  system  is  normal  to  the  local  velocity  vector 
because  an  LV  with  only  one  set  of  fringes  responds  to  only  one  component 
of  the  velocity  at  a  point.  Since  there  are  very  few  LV’s  that 
simultaneously  respond  to  all  three  possible  components  of  the  velocity 
vector  and  turbulent  flows  are  three  dimensional  in  the  fluctuation 
components  of  velocity,  a  more  general  relationship  for  turbulent  flow 
measurements  is  useful.  Consider  a  turbulent  flow  field  in  which  the 
mean  motion  can  be  characterized  primarily  by  a  single  component,  U.  The 
velocity  vector  at  a  point  may  then  be  written  as 

Vj  =  (<U>  +  uf)  i  vj  +  wfk  (24) 


where  u  ,  v , ,  w  ,  are  the  fluctuation  components  along  the  coordinate 

"h  ~t  -+ 

axes  defined  by  the  orthogonal  unit  vectors  (r,  3,  k) .  The  magnitude 
of  the  velocity  vector  is  then 


VJ  =  <U  >  [  1 


<U> 


u2  +  v2  +  w2  \2 

-^>2  ') 


(25) 


Substituting  Eq.  (25)  in 
axial  component  measured 


Eq.  (15)  with  f  =  <U>  + 
by  an  LV  system 


(1  +  ^-) 

\  <u>/ 


1=1 


U  =  <U> 


(l+O,) 


M 

£— L- 
1=  1 

<i  +  a:> 


gives  for  the  mean 


(26) 
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where 


=  2 


u2  + 


2  2 
v  +  w 


<u> 


<  U  >2 


(27) 


Expanding  Eq.  (26)  in  a  binomial  series  and  truncating  terms  of  third  or 
higher  order  leads  to 


<U  > 


U  /1+g<'>  ~  1  (g<v>  + 


(28) 


Figure  3  predicts  the  bias  expected  in  turbulent  flow  IRLV  meas¬ 
urements  when  u  -a  ~  2/ 3a  which  Sami  (Ref.  38)  indicates  are 
<v>  <w>  <u> 

representative  of  the  near  field  of  a  subsonic  jet.  The  bias  predicted 


Figure  3.  Bias  of  the  axial  velocity  component  in  turbulent  flows  (S  -  0). 
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for  the  arithmetically  averaged  component  velocity  is  obviously  less 
than  the  bias  that  would  result  in  a  one-dimensional  flow.  But,  it 
is  still  large  enough  to  be  significant  in  many  IRLV  applications.  Note 
that  McLaughlin  and  Tiederman  (Ref.  19)  and  Buchhave  (Ref.  24)  also  show 
component  bias  to  be  generally  less  in  three-dimensional  than  in  one¬ 
dimensional  flow. 

Note  that  Eq.  (28)  is  similar  to  the  expression  derived  for  a  purely 
Dne-dimensional  flow,  Eq.  (20),  when  the  turbulent  intensities  associated 
with  the  (j ,k)  axes  become  negligible.  However,  the  turbulent  intensity 
associated  with  the  bias  in  Eq.  (28)  is  the  arithmetic  value  rather  than 
the  temporal  turbulent  intensity.  This  is  a  consequence  of  the  averaging 
techniques  and  truncations  used  in  the  derivation;  but,  as  noted  from  Fig. 

3,  the  error  involved  in  the  estimate  is  not  large  for  turbulence  levels 
of  30  percent  or  less.  Since  the  binomial  series  will  only  converge  for 
small  values  of  o^,  taking  a<u>  =  a-  in  Eq.  (28)  is  not  a  restrictive 
assumption. 

2.2  EXPERIMENTAL  SIMULATION 

With  a  relationship  for  the  bias  magnitude  derived,  the  flow  conditions 
desired  in  an  experimental  study  of  bias  can  be  delineated.  Specifically, 
the  experiment  must  have  (1)  constant  fluid  density,  (2)  spatially  and 
temporally  invariant  particle  number  density,  (3)  negligible  particle- 
fluid  lags,  (4)  single  counting  of  each  particle  (i.e.,  rapid  data 
acquisition  system  and  negligible  fluid  recirculation)  and  (5)  particles 

A 

entering  the  probe  volume  frequently  enough  that  is  approximately 


An  air  jet  was  selected  for  the  experiments  to  provide  a  wide  range 
of  turbulent  intensities  in  all  component  directions.  Thus,  arithmetic 
means  of  IRLV  one-component  data  should  be  sufficiently  biased  for 
experimental  identification  in  at  least  part  of  the  flow.  The  jet  exit 
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Mach  number  was  kept  below  0.1  to  ensure  that  the  fluid  density  could 
be  considered  constant.  A  fairly  constant  particle  number  density  was 
obtained  by  filtering  the  supply  air  to  remove  all  detectable  particles 
and  then  seeding  with  a  fluidized  bed  seeder  (Fig.  4),  Additionally, 
the  jet  was  enclosed  so  that  recirculation  external  to  the  jet  tended 
to  spatially  equalise  particle  concentration. 

By  varying  the  supply  pressure  to  the  fluidized  bed  seeder  the  seed 
concentration  was  controlled  and  a  range  of  particle  interarrival  times, 
Ati#  obtained.  Based  on  particle  dynamic  calculations  for  the  range  of 
conditions  in  the  jet,  uniform  1— pm-diam  particles  were  used  in  the 
seeder  because  particles  of  that  size  accurately  respond  to  fluid 
transients.  Fluidized  bed  seeders  will  not  provide  a  perfectly  uniform 
particle  concentration  (Ref.  39);  however,  the  vibrating  seeder  shown  in 
Fig.  4  provided  a  more  uniform  particle  number  density  than  conventional 
seeders.  The  jet  was  also  seeded  with  a  Berglund-Liu  raonodisperse 
aerosol  generator  (Ref,  40)  in  a  few  of  the  tests,  but  this  seeder 
could  produce  only  low  seed  concentrations.  Air  from  the  seeder  entered 
the  jet  stilling  chamber  after  passing  through  a  particle  charge  neu¬ 
tralizer.  The  neutralizer  lessened  particle  agglomeration  by  min¬ 
imizing  electrostatic  attraction,  but  some  coagulation  was  observed  in 
photomicrographs  of  particles  collected  from  the  test  flow. 

The  stilling  chamber  and  jet  test  apparatus  is  shown  in  Fig.  5.  A 
coaxial  mixer  section  of  the  stilling  chamber  mixed  the  seeded  air  with 
filtered  air  to  establish  the  desired  seed  concentration.  A  shielded, 
self-aspirating,  copper-constantan  thermocouple  and  two  static  pressure 
ports  in  the  center  of  the  stilling  chamber  monitored  the  flow  total 
conditions.  A  25-to-1  area  ratio  contraction  nozzle  with  a  2.54-cm- 
diam  throat  was  at  the  end  of  the  stilling  chamber.  The  nozzle  con¬ 
tour  was  matched  to  calculated  potential  flow  streamlines  to  smoothly 
accelerate  the  air  and  particles  to  a  uniform  nozzle  exit  velocity. 
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The  jet  was  enclosed  in  a  recirculation  test  chamber  that  was  18 
nozzle  diameters  long  and  5  diameters  square.  The  front  wall  of  this 
chamber  was  a  thin,  optically  flat  window  which  allowed  LV  measurements 
at  axial  stations  up  to  X/D  =  9.  The  other  chamber  walls  were  painted  a 
flat  black  to  minimize  light  reflections  that  could  interfere  with  LV 
signals.  Static  pressure  ports  along  the  back  wall  of  the  test  chamber 
were  at  stations  X/D  =  2,  4,  and  8.  A  slot  in  the  bottom  of  the 
chamber  could  be  opened  to  insert  pitot  and  hot-wire  probes  mounted  on 
the  two-component  traverse  mechanism  below  the  chamber. 

2.2.1  LV  Instrumentation 

The  LV  system  and  associated  electronics  used  in  this  study  have 
been  described  in  detail  in  Refs.  1  and  2.  The  present  description 
will,  therefore,  concentrate  on  those  factors  influencing  overall 
LV  accuracy. 

To  measure  statistical  bias  experimentally,  the  precision  uncertainty 
in  the  IRLV  measurements  must  be  less  than  the  bias.  The  IRLV  and  traverse 
system  shown  in  Fig.  6  provided  this  precision  in  the  low-speed  jet 
measurements.  A  514.5-nm  wavelength  light  beam  from  an  argon- ion  laser 
was  diffracted  into  several  beams  by  the  water-filled  Bragg  cell.  The 
Bragg  cell  beam  splitter  described  in  Refs.  1,  2,  and  15  measured  neg¬ 
ative  velocities  which  frequently  occur  in  highly  turbulent  jet  flows. 

One  nearly  equal  intensity  beam  pair  from  the  Bragg  cell  was  passed  through 
the  collimating  lens,  the  beam  expander,  and  the  redirecting  optics 
before  being  focused  to  form  the  LV  probe  volume.  Although  this  system 
could  measure  two  velocity  components  (see  Ref.  2),  only  one  was  measured 
to  simplify  and  speed  up  data  acquisition. 

For  a  dual-scatter  LV  system,  Refs.  1  and  6  give  the  precision 
uncertainty  as 

in,  -  a  -  $u(  +  yu?  (29) 
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Figure  6.  Laser  velocimeter  and  traverse  system. 
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where 

«'-  K8f„  (*.  .  I.A,.) 

A  Kt 

B  =  _ 5.  -  2f  At 

“  T-  ok 


K.  is  the  LV  fringe  spacing,  f  is  the  Bragg  system  reference  frequency, 

6  o 

and  t  is  the  measured  LV  burst  period.  The  uncertainty  in  the  LV  fringe 

spacing,  AK^/K^,  evaluated  as  described  in  Ref.  1  from  LV  calibrations, 

is  +0.037  for  the  system  K„  of  18.3  pm.  The  uncertainty  in  the  1 5-MHz 

A  ^ 

reference  frequency  is  estimated  as  AfQ/f  =  +6.66  x  10  ,  and  the 
uncertainty  of  the  LV  burst  period,  was  determined  in  Ref,  2  to  be 

+0.005  nsec  for  the  particular  LV  burst  counter  processor  used. 

Resulting  precision  uncertainties  for  velocities  measured  by  this  LV 
system  are  summarized  in  Table  1 . 

Table  1.  Precision  Uncertainty  of  Laser  Velocimeter 


(30) 

(31) 

(32) 


Velocity,  U 

Precision,  AU 

m/sec 

ft/sec 

ra/sec 

ft/sec 

0 

0 

±0.20 

±0.68 

±  7.6 

±  25 

±0.48 

±1.58 

±15.6 

±  50 

±0.76 

±2.50 

±22.9 

±  75 

±1.05 

±3.43 

±30.5 

±100 

±1.33 

±4.35 

±38.1 

±125 

±1.61  ] 

±5.28 
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Because  a  finite  number  of  measurements  determines  IRLV  averages, 
the  statistical  confidence  of  these  measurements  must  be  considered. 
The  95-percent  confidence  interval  for  the  arithmetic  mean  velocity  is 
given  in  Ref.  1  as 


A  <  V  > 


(33) 


The  95-percent  confidence  interval  for  the  turbulence  intensity  is 
given  in  Ref.  41  as 


A<y 


<i'> 


=  ml 


.282 

\  M  4 


(34) 


where  K  is  the  kurtosis  coefficient,  defined  as 

K  =  L  f  Tv^V^ 


dt 


-I 

2 


(V  a*) 


(35) 


To  establish  a  better  than  1.5-percent  confidence  interval  for  the  mean 
velocity,  all  averages  were  calculated  from  at  least  1,500  samples.  For 
this  sample  size,  a  better  than  3-percent  confidence  interval  resulted 
for  the  turbulence  intensities. 


The  uncertainty  in  positioning  the  IRLV  probe  volume  will  also 
contribute  to  uncertainty  in  the  IRLV  measurements.  With  the  traverse 
system  shown  in  Fig.  6  the  probe  volume  could  be  moved  along  any  of 
three  mutually  orthogonal  axes  and  positioned  to  within  +0.08  cm  (+0.03 
in.)  by  properly  accounting  for  hysteresis.  However  there  is  an  ad¬ 
ditional  uncertainty  in  locating  the  LV  measurements  due  to  the  finite 
size  of  the  probe  volume.  The  diameter  of  the  probe  volume  is  taken  as 
the  diffraction  limit  of  the  focusing  lens.  For  the  0.381-m  focal  length 
lens  and  1 . 5  mm  beam  diameters  at  the  lens  a  probe  volume  diameter  of 
0.17  mm  results.  The  apparent  probe  volume  length  was  reduced  to  ap¬ 
proximately  0.33  cm  by  viewing  the  probe  volume  through  a  one-mm-di- 
ameter  pinhole  at  20  deg  off-axis  in  backscatter. 
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The  LV  signal  processing  and  recording  systems  achieved  the  accuracy 
required  to  measure  the  anticipated  statistical  bias,  while  attaining 
near-maximum  IRLV  data  rates  (and  minimum  <At>)  for  this  experiment. 
Signals  from  the  photomultiplier  (PMT)  were  processed  without 
amplification  to  eliminate  processing  the  weak,  noisy  signals  from  the 
edges  of  the  probe  volume.  Also  the  PMT  signals  were  processed  only 
when  the  pedestal  of  the  burst  was  detected  by  the  high-frequency  re¬ 
jection  A-gate  oscilloscope  trigger.  Then  the  fringe-crossing  period 
of  each  triggered  signal  was  determined  by  a  Model  8  Doppler  data 
processor  (DDP)  only  if  the  signal  passed  the  processor  noise  rejection 
criteria.  These  criteria  required  that  the  fringe— crossing  periods 
based  on  4,  5,  and  8  cycles  agree  to  within  1.5  percent,  and  only  signals 
between  8  and  20  MHz  (-76  to  -122  m/sec)  were  processed  to  reduce  the 
possibility  of  recording  noise. 

Additional  details  of  the  model  8  DDP  are  given  in  Refs.  15  and  31 . 
The  processor  used  was  modified  by  interval-counting  circuitry  triggered 
by  the  DDP  print  command  to  measure  the  time,  At^,  between  successive 
validated  signals.  The  interval-counting  resolution  was  +0 . 5  psec. 

With  this  modification,  approximate  time-integrated  average  velocity 
moments  and  the  data  rates  in  addition  to  the  arithmetic  average  moments 
were  determined. 

To  attain  near-maximum  IRLV  data  rates  the  data  were  recorded  using 
direct  memory  access  to  a  minicomputer.  This  system  could  typically 
record  over  13,000  fringe  periods  and  time  intervals  per  second.  At  this 
rate  and  a  15— m/sec  mean  velocity  there  was  a  high  probability  that  the 
LV  signals  would  result  from  mutliple  particle  traversals  of  the  probe 
volume,  but  might  not  be  processed  by  the  DDP.  By  controlling  the  seeder 
pressure,  the  flow  was  never  seeded  heavily  enough  to  exceed  mean  data 
rates  of  4,000  per  second;  this  lessened  the  possibility  of  excluding 
measurements  because  of  multiple  particle  traversals. 
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2.2.2  Supporting  Instrumentation 

Statistical  biasing  will  not  occur  in  averages  of  velocities  measured 
by  continuously  monitoring  the  jet  flow.  Such  continuous  measurements 
were  made  using  both  pressure  and  hot-wire  instrumentation  and 
statistical  bias  was  identified  by  comparing  the  average  LV  measurements 
with  the  average  pressure  and  hot-wire  measurements.  However,  such 
identification  is  possible  only  if  the  uncertainty  in  the  pressure  and 
hot-wire  measurements  is  less  than  the  bias.  Furthermore,  the  pressure 
and  temperature  stilling-chamber  measurements  used  to  monitor  the  jet 
flow  must  be  equally  accurate  for  a  meaningful  comparison  of  the  in¬ 
dependent  measurements. 

To  attain  the  necessary  accuracy  in  continuous  pressure  measurements 
Barocel  pressure  transducers  were  used.  These  transducers  maintained  an 
accuracy  of  +0.15  percent  of  the  measured  pressure  throughout  their 
+1 0—mm-Hg  pressure  range,  as  determined  by  calibration  against  a 
laboratory  standard  Hook  gage  which  was  traceable  to  the  National 

Bureau  of  Standards  (NBS) .  The  ambient  pressure,  P  ,  was  measured  with 

3  • 

a  Bourdon-tube  pressure  gage,  and  the  air  temperature,  T,  was  monitored 
using  a  shielded,  self-aspirated,  copper-constantan  thermocouple.  Both 
instruments  were  calibrated  against  NBS  laboratory  standards.  The 
measurements  were  continuously  monitored  by  the  minicomputer  data  system 
which  continuously  displayed  the  value  of  each  variable  as  well  as  the 
velocity  calculated  from 


where  g  and  J  are  unit  conversion  factors,  c  is  the  fluid  specific  heat. 
•  c  c  p  7 

T  is  the  jet  fluid  temperature,  P  is  the  ambient  pressure,  P  is  the  test 

a  c 

cell  gage  pressure,  AP  is  the  stilling  chamber  or  pitot  probe  to  test  cell 
differential  pressure,  and  y  is  the  ratio  of  specific  heats.  The 
precision  of  these  velocity  measurements  was  determined  by  expanding  Eq. 
(36)  in  a  Taylor  series  in  all  of  the  measured  variables  and  evaluating 
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the  expression  at  the  maximum  uncertainty  of  each  (see  Ref,  42) .  The 
analysis  yielded  an  uncertainty  in  velocity  of  less  than  1  .5  percent  at 
the  nominal  jet  exit  velocity  of  30.5  m/sec. 

The  hot-wire  measurements  were  made  using  a  Thermo-Systems,  Inc. 
1050  series,  constant-temperature  anemometer  equipped  with  a  3.8-pm- 
diam,  tungsten  sensor.  The  ends  of  the  1.5— mm-long  sensor  were  plated 
to  reduce  end  effects,  leaving  an  active  sensor  length  of  1.25  mm. 
Assuming  the  sensor  would  properly  respond  only  to  eddies  of  this  size 
or  larger,  the  frequency  response  of  the  sensor  was  estimated  to  be 
24  kHz  when  eddies  were  convected  past  the  sensor  at  the  nominal  jet 
exit  velocity.  Therefore,  signal  frequency  components  greater  than 
30  kHz  were  considered  noise  and  filtered  out.  The  hot-wire  bridge 
output  was  recorded  on  FM-modulated  analog  tape  maintaining  a  30-kHz 
frequency  response. 


To  attain  the  needed  accuracy  in  the  hot-wire  measurements,  the 
recorded  bridge  output  was  digitized  at  a  50-kHz  rate  and  each  sample 
was  converted  to  a  velocity  using  the  hot-wire  response  equation  (Refs. 
43  and  44), 


e?  =  A  -  Bl'|  CU, 


(37) 


The  constants  A,  B,  C,  and  h  were  determined  by  a  nonlinear  least— squares 
fit  of  Eq.  (37)  to  the  hot-wire  calibration  (Ref.  45).  Calibrations 
were  performed  every  30  minutes  in  the  test  apparatus  using  velocities 
measured  with  the  pressure  instrumentation.  The  calibration  data 
exhibited  excellent  agreement  during  all  but  one  test,  and  those  test 
data  were  not  used  for  the  comparisons  shown.  Figure  7  shows  two 
typical  calibrations.  Because  Eq.  (37)  fit  the  calibration  data  so  well 
(Fig.  7)  the  uncertainties  in  the  hot-wire  velocity  measurements  were 
only  slightly  larger  than  the  uncertainties  in  the  pressure-determined 
calibration  velocities.  Furthermore,  by  using  approximately  200,000 
digital  hot-wire  samples  when  averaging  the  hot-wire  data,  the  95-percent 
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statistical  confidence  intervals  were  better  than  0.1  percent  of  the  mean 
velocities  and  approximately  0.2  percent  of  the  turbulence  intensities 
and  should  not  add  significantly  to  the  hot-wire  measurement  uncertainty. 
The  hot-wire  and  pitot  probe  position  uncertainties  were  +0.5  mm. 

2.3  RESULTS 

The  experimental  examination  of  statistical  bias  was  conducted 
in  three  phases.  First,  the  jet  flow  was  evaluated  to  determine  how 
well  the  flow  conditions  that  theoretically  result  in  bias  were  met. 
Second,  the  flow  was  mapped  with  pitot  and  hot-wire  probes.  And,  third, 
the  flow  was  mapped  with  the  IRLV.  The  second  and  third  phase  results 
were  then  compared  to  determine  the  magnitude  of  any  bias  present. 


Figure  7.  Hot-wire  anemometer  calibration. 


34 


AEDC-TR  -78-65 


In  the  first  phase  the  operation  of  the  seeder  was  examined  to 
determine  whether  it  provided  a  steady,  regular,  controllable  output. 
Successive  LV  detections  of  particles  in  the  jet  core  were  recorded, 
and  the  time  between  each  detection,  At^,  was  measured  by  the  LV  data 
processor.  Histograms  of  these  At^  data,  presented  in  Fig.  8  for 
several  seeder  settings,  show  that  nearly  equal  times  separate  most 
particle  detections;  the  standard  deviation  of  the  times  between 
detections  is  less  than  0.5  percent  of  the  average  time  between  detections 
for  all  cases.  In  the  smooth  core  flow  these  nearly  equal  times  indicate 
nearly  equal  fluid  mass  separating  each  particle,  and  thus  a  fairly 
uniform  particle  concentration.  Note  that  no  particles  were  detected  with 
the  seeder  turned  off,  verifying  that  only  particles  from  the  seeder 
comprise  the  LV  data  sets  obtained  in  the  tests. 


The  distribution  of  the  seed  concentration  in  the  test  chamber  was 

also  examined  by  timing  successive  LV  detections  and  evaluating  the 

average  time  between  detections,  <At>.  If  the  ratio  <At  >/<At(r)>  is 
-  E 

distributed  as  U(r)/U„  throughout  the  jet,  then  it  can  be  concluded 

from  Eq.  (9)  that  the  particle  concentration  is,  on  the  average,  uniform 

throughout  the  jet.  Figure  9  shows  that  radial  distributions  of  <At„>/ 

_  E 

<At(r)>  do  correspond  to  pitot  probe  measurements  of  u(r)/u„,  illustrating 

E 

that  the  particle  concentration  is  essentially  uniform,  on  the  average, 
throughout  the  jet. 


To  examine  the  size  of  particles  in  the  test  flow,  particle  samples 
were  collected  near  the  test  chamber  exit.  Photomicrographs  taken  of 
these  samples  using  a  scanning  electron  microscope  show  that  some  coagulation 
of  the  1-ym  seed  occurred.  The  distribution  of  particle  sizes  determined 
from  several  photomicrographs  is  shown  in  Fig.  10  and  indicates  that  most 
particles  had  1-  to  3-ym  diameters.  A  few  larger  particles  were  visable 
in  the  photomicrographs,  but  they  were  chain-shaped  coagulations  that  most 
likely  were  created  during  the  particle  collection.  Particle  dynamic 
calculations  indicated  that  3-ym  particles  did  not  appreciably  lag  the 
flow,  so  the  departure  from  the  desired  1-ym  size  had  little  effect 
on  the  results. 
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Probability  P(At) 


a.  Average  rate,  25.8  samples/sec  b.  Average  rate,  343.9  samples/sec  c.  Average  rate,  3235.8  samples/sec 

(note  scale  change) 


Figure  8.  Histograms  of  particle  interarrival  times. 


AEDC-TR  -78-65 


Average  Rate,  <fitr  > /<&t{r) 


a.  Variation  of  data  rate  at  X/D  =  4 
Figure  9.  Data  rate  variation  in  jet  near  field. 
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Average  Rate,  <At->/<at(r) 
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Measured  Particle  Diameter,  Dp,  pm 

Figure  10.  Particle  size  distribution. 


Once  the  seeding  characteristics  were  determined,  the  flow  field 
was  mapped  by  conventional  instrumentation.  Radial  surveys  were  made 
at  axial  stations  of  2,  4,  6,  and  S  nozzle  diameters  with  both  pitot 
and  hot-wire  probes.  A  hot-wire  radial  survey  at  X/D  *  1/16  measured 
the  jet  core  turbulence,  which  was  less  that  2  percent  of  the  nominal 
exit  velocity.  Figure  11  shows  the  results  of  the  surveys  of  the  mean 
axial  velocity  where  the  mean  velocities  were  fit  with  the  semi- 
empirical  expression 


U  iX,  rl 


exp 


-0.693 


VV)  V 

V  MX)  / 


(38) 


as  an  aid  for  comparison  to  the  LV  results. 
U (X) ,  inner  radius  of  the  jet  shear  layer. 


The  centerline  velocity, 
r  (X) ,  and  the  shear  layer 


width  to  the  half-velocity  point,  b(X) ,  were  determined  using  the  curve¬ 
fitting  method  of  Ref.  45  and  are  given  in  Table  2.  Even  though  only  hot- 
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wire  data  points  are  indicated  in  Fig.  11,  both  hot-wire  and  pitot 
measurements  were  used  in  deriving  the  curve  fits. 

Although  measurements  were  made  across  the  entire  jet,  only  data 
from  the  jet  region,  where  the  local  turbulence  intensity  was  less  than 
0.33,  will  be  used  for  bias  identification.  Data  aliasing  attributable 
to  negative  velocities  and  hot-wire  response  variations  can  be  sig¬ 
nificant  in  pitot  and  hot-wire  measurements  at  higher  intensities  (Ref. 

37  and  46).  Note  that  hot-wire  measurements  were  only  made  in  unseeded 
flows.  The  hot-wire  calibration  tended  to  drift  more  in  seeded  flows, 
and  particle— wire  impingements  can  cause  pseudo-velocity  signals  that 
alter  the  hot-wire  accuracy.  However,  pitot  measurements  made  in  the 
seeded  flow  verify  that  the  mean  velocity  characteristics  of  the  jet  at 
high  particle  concentrations  were  identical  to  those  of  the  unseeded  flow. 

The  turbulence  intensities  measured  with  the  hot  wire  at  axial  stations 
of  2,  4,  6,  and  8  nozzle  diameters  are  shown  in  Fig.  12.  The  data  inside 
the  a-  =  0.33  limits  were  least  squares  fitted  to  a  second  order  polynomial: 

l  A  (x)  +  B  (x)  r  +  C(x)r2  ,  for  r  <  r 
_  J  o  c t  _  a 

}  A  (x)  +  B(x)r  +■  C(x)i2  ,  for  r  >  j  (39) 

The  polynomial  parameters  A(x) ,  B(x),  C(x),  and  r  are  given  in  Table  2. 

CT 

The  polynomial  data  fits,  along  with  Eqs.  (20)  and  (38),  predicted  the 
maximum  mean  velocity  bias  that  could  be  anticipated  in  the  IRLV  meas¬ 
urements.  The  anticipated  bias  is  shown  in  Fig.  11  along  with  the  hot¬ 
wire  uncertainty.  In  determining  the  hot-wire  uncertainty,  the  hot¬ 
wire  pressure  calibration  uncertainty,  derived  from  Eq.  (36),  was  doubled 
to  account  for  any  other  hot-wire  measurement  uncertainties.  Even  with 
the  broad  uncertainty  assumption  it  is  obvious  from  Fig.  11  that  the  IRLV 
statistical  bias  should  be  identifiable  in  these  tests. 


40 


Table  2.  Jet  Flow-Field  Data 


X/D 

~i 

Mean  Velocity  Fit 

Eq.  (38) 

Turbulence  Intensity  Fit  Eq.  (39) 

(x) 

rQ(x) 

b  (x) 

A(x) 

B(x) 

C  (x) 

ra 

0.06 

99.14 

0.479 

0.021 

2.0 

98.97 

0.235 

0.319 

0.0558 

-0.3976 

1.7308 

1.1831 

4.0 

96.07 

0 

0.563 

0,0896 

0.1644 

0.6238 

0 

6.0 

83.57 

0 

0.632 

0.1606 

0.2072 

0.2865 

0 

o 

a 

00 

66.98 

0 

0.737 

0.2332 

0.1418 

0.2385 

0 

Velocity,  u/u 


Radius,  r/D 


a.  X/D  =  1/16 

Figure  11.  Hot-wire  measurements  of  axial  velocity  distributions. 
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Velocity,  u/u 


b.  X/D  =  2 
Figure  11.  Continued. 
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Velocity,  u/u 


c.  X/D  =  4 

Figure  11.  Continued. 


AEDC-TR-78-65 


Velocity,  u/u 


d.  X/D  =  6 
Figure  11.  Continued. 
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Velocity,  u/u 


e.  X/D  =  8 
Figure  11.  Concluded. 
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In  the  theoretical  discussion  of  statistical  bias,  requiring  -  V± 
was  equivalent  to  assuming  that  an  IRLV  would  detect  several  particles 
per  integral  time  scale;  bias  might  be  reduced  if  detections  rates  are 
not  this  high.  To  examine  this  assumption,  integral  time  scales  were 
determined  from  the  digital  hot-wire  data  using  a  fast  Fourier  transform 
technique  (Ref.  47).  To  reduce  computer  storage  requirements,  hot-wire 
measurements  acquired  over  a  15-second  interval  were  segmented  into  small 
digital  data  records  of  about  0.17  second  duration  and  transformed 
record  by  record.  The  resulting  power  spectra  wer£_averaged  and  trans¬ 
formed  to-determine  the  autocorrelations  from  which  the  integral  time 
scales  were  determined  by  numerical  integration.  Integral  time  scales 
obtained  in  this  manner  are  compared  to  values  interpolated  from  the 
data  of  Sami  (Ref.  38)  in  Fig.  13.  Data  rates  of  1,500  to  5,000  samples  per 
sec  were  necessary  to  ensure  one  particle  per  Integral  length  scale  at 
axial  stations  between  X/D  =  4  and  8  where  the  bias  was  to  be  identified. 
Thus,  to  test  for  bias-data-rate  dependence,  data  were  recorded  at  rates 
ranging  from  10  to  4,800  samples/sec. 

In  addition,  the  hot-wire  data  were  used  to  determine  the  radial 
variation  of  the  skew  and  the  kurtosis  shown  in  Figs.  14  and  15,  respec¬ 
tively.  These  higher  order  central  moments  were  used  for  predicting  the 
statistical  bias  in  IRLV  averages  of  the  turbulence  intensity  and  skew 
coefficient  as  indicated  by  Eq.  (22). 

After  the  "unbiased"  flow-field  characteristics  were  established, 

IRLV  radial  surveys  were  made  at  axial  stations  of  2,  4,  6,  and  8  nozzle 
diameters.  Identification  of  statistical  bias  primarily  involved 
examination  of  the  X/D  =  4  and  6  data.  At  X/D  =  2  the  shear  layer  had 
enough  curvature  to  produce  mean  velocity  gradients  in  the  LV  probe  volume, 
making  those  measurements  questionable.  At  X/D  =  8  the  mean  velocity 
was  low,  and  the  pitot  and  hot-wire  accuracies  were  not  sufficient  for 
bias  Identification  throughout  the  flow. 
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Radial  Location,  r/D 

Figure  13.  Radial  variation  of  integral  time  scales. 

The  IRLV  radial  surveys  were  repeated  at  several  seeder  settings 
to  define  any  bias  variation  with  seed  density.  These  mean  velocity 
surveys  are  compared  in  Fig.  16  with  the  "unbiased"  velocity  curves. 
Note  that  the  IRLV  data  more  closely  approximate  the  unbiased  velocity 
curves  than  the  theoretical  bias  curves,  regardless  of  the  seeding  rate. 
The  IRLV  measurements  of  turbulence  intensity  (Fig.  17)  also  agree  well 
with  the  unbiased  velocity  curves,  indicating  that  statistical  biasing 
was  not  nearly  as  significant  as  predicted  for  these  tests.  No  sig¬ 
nificant  bias  dependence  on  data  rate  was  discerned. 

To  date  the  only  AEDC  tests  that  show  differences  in  IRLV  and 
hot-wire  measurements  that  might  be  attributable  to  bias  were  reported 
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Skew  Coefficient 


Figure  14.  Radial  distribution  of  skew  coefficient. 
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Kurtosis 


Figure  15.  Radial  distribution  of  axial  kurtosis  coefficient. 
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Figure  16.  Comparison  of  IRLV  mean  velocity  with  "biased"  predictions. 
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Figure  17.  Comparison  of  IRLV  and  hot-wire  local  turbulence  intensities. 
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in  Ref.  11.  In  that  experiment,  the  difference  in  the  two  measurements 
was  most  obvious  in  the  centerline  axial  velocity  decay.  Accordingly, 
the  centerline  axial  velocity  decay  and  turbulence  intensity  distribution 
as  determined  by  pitot,  hot-wire,  and  IRLV  data  averages  in  the  present 
experiment  are  plotted  in  Fig.  18.  These  data  show  no  evidence  of  either 
statistical  bias  or  a  bias-data-rate  dependence. 
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Figure  18.  Centerline  distribution  of  axial  velocities 
and  turbulence  intensities. 

In  this  study  bias  was  identified  by  comparing  averages  of  IRLV 
data  with  unbiased  measurements.  It  is  clear  from  Eq.  (6)  and  stated 
in  Ref.  20  that  if  the  particle  interarrival  times  are  random  no  bias 
will  result.  Reference  25  suggests  further  that  random  selection  of 
samples  from  an  IRLV  data  set  will  lead  to  unbiased  averages  of  the 
velocity  moments.  Therefore,  in  a  final  attempt  to  identify  statistical 
bias,  the  IRLV  averaged  mean  velocities  and  turbulence  intensities  were 
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compared  with  averages  of  randomly  selected  data.  The  selection  was 
done  several  times,  first  selecting  one  sample  at  random  from  each  200 
recorded  samples,  then  one  from  each  600,  and  finally  one  from  each 
2,000  recorded  samples.  The  ratio  of  the  original  IRLV  data  average  to 
the  randomly  selected  data  average  for  various  radial  locations  at  X/D 
of  4  and  6  is  shown  in  Fig.  19a.  The  ratio  of  the  original  IRLV  data 
root-mean-square  average  to  the  randomly  selected  root-mean-square  for 
the  same  X/D  locations  is  shown  in  Fig.  19b.  There  is  no  evidence  of 
statistical  bias  in  either  data  set. 

2.4  ANALYSIS  OF  INFLUENTIAL  FACTORS 

Based  on  the  preceding  results,  the  magnitude  of  IRLV  statistical 
bias  can  be  much  less  than  predicted  by  the  theory  outlined  previously. 
Several  factors,  which  could  have  reduced  the  bias  magnitude  in  the 
preceding  results,  will  be  examined  to  determine  if  the  lack  of  bias 
might  be  unique  to  these  results.  The  factors  examined  include  (1) 
the  relationship  between  sampling  rate  and  turbulence  structure,  (2) 
processing  rates  possible  with  a  given  data  acquisition  system,  and 
(3)  fluctuations  in  the  local  particle  concentration. 

2.4.1  Sampling  Rate 

In  Ref.  20,  Barnett  and  Bentley  concluded  that  the  bias  between 
arithmetic  and  time-integrated  averages  of1  IRLV  data  will  decrease  as 
the  sampling  rate  per  primary  flow  frequency  goes  down.  This  conclusion 
can  be  supported  using  the  previously  presented  theory  by  evaluating  the 
case  when  the  time  interval  between  detections,  At^,  is  large  compared  to 
the  flow  integral  time  scale,  LT.  In  this  limiting  case,  as  the  average 
interval  time  to  integral  scale  ratio,  <fit>/L  ,  becomes  large,  the 
velocity  average  between  detections,  V^,  will  approach  the  temporal  mean 
velocity,  V,  and,  from  Eq.  (3)  arithmetic  averages  will  approach  ap¬ 
proximate  temporal  averages.  As  a  result,  if  the  approximate  temporal 
averages  are  still  good  approximations  of  the  true  temporal  averages  as 
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Data  Set  Mean  Ratioed  to  Randomly  Sampled  Set  Mean 


a.  Ratio  of  velocity  means 
Figure  19.  Effect  of  data  randomization. 
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Means  Squares  Ratio  (Randomly  Sampled/Original  Data) 


b.  Ratio  of  velocity  variance 
Figure  19.  Concluded. 
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the  time  between  detections  becomes  much  larger  than  the  integral  time, 

L^,  then  statistical  biasing  must  eventually  decrease  with  increasing 
time  between  detections  and  finally  disappear. 

To  examine  the  possible  dependence  of  statistical  bias  on  data  rate, 
IRLV  signals  were  simulated  in  a  numerical  experiment.  In  this  simulation 
a  function  was  generated  to  represent  the  flow  velocity  in  the  probe 
volume  at  time  t.  Velocities  were  then  recorded  as  particle  realizations 
every  time  a  specified  amount  of  fluid  entered  the  probe  volume.  For 
consistency  with  the  bias  theory,  negative  velocities  were  not  allowed  to 
recirculate  particles  through  the  probe  volume.  Initially,  a  sinusoidal 
velocity  V  =  Vq  +  sin  (tot  +  q>)  was  assumed  for  comparison  with  the 
part icle-per-main-f low  frequency  criteria  of  Ref.  20.  The  quantity  At/L^, 
given  in  this  case  as  At/(1/w),  was  varied  by  changing  either  ui  or  the 
specified  fluid  mass  separating  particles.  Both  the  arithmetic  average 
velocity  and  the  approximate-time-integ rated  average  velocity  were  compared 
with  the  true  time-averaged  velocity. 

The  results  of  the  simulation  are  shown  in  Fig,  20.  Since  all 

theorized  statistical  bias  conditions  are  met  perfectly  in  this  simulation, 
*  2 

the  quantity  (<V>/V  -  1)/cr^  should  vary  between  unity  for  completely 
biased  cases  and  zero  for  unbiased  cases.  Excursions  beyond  these  limits 
result  when  the  sampling  occurs  as  a  multiple  of  the  flow  frequency. 

As  indicated  in  Fig.  20a,  the  bias  in  arithmetic  averages  is  not 
reduced  when  the  average  time  between  detections  is  large  (i.e. ,  the 
data  rate  per  flow  frequency  is  small)  as  suggested  above.  However,  in 
the  above  analysis  bias  will  be  eliminated  only  when  the  time  between 

2Lt 

detections  is  large  enough  to  drive  the  product  (V^-V)  in  Eq.  (12) 

toward  zero.  For  the  data  shown  in  Fig,  20a,  the  lowest  value  of  the 

2  /  ^^^^****2.  2  2 

product  — ; (V.-V)  is  greater  than  18  ft  /sec  .  To  decrease  this  prod- 

<ut>  1 

uct  to  near  zero  (three  orders  of  magnitude)  requires  a  corresponding 
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Fraction  of  Estimated  Bias, 


a.  Arithmetically  averaged  data,  <v> 

Figure  20.  Simulation  program  predictions  of  sampling  rate  effect  on  the  velocity  bias  of  a  sinusoidal  flow. 
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decrease  in  the  data  rate  to  about  20  samples  per  hour.  In  the  simulation 
the  bias  quantity  (<V>/V  -  1)/o*  is  reduced  to  only  about  0.8  at  data 
rates  approaching  20  samples  per  hour,  and  even  lower  data  rates  are 
required  for  completely  unbiased  arithmetic  averages.  These  lower  data 
rates  will  be  much  too  low  for  practical  use  of  a  laser  velocimeter. 
Simulation  of  the  impractical  lower  data  rates  is  not  expedient  because 
the  large  values  of  the  average  time  interval  require  correspondingly 
long  computer  run  times. 

Figure  20b  shows  that  at  values  less  than  one  sample  per  flow 
frequency  the  approximate  time-integrated  velocity  is  biased.  This 
verifies  Barnett  and  Bentley's  hypothesis  that  arithmetic  averaging 
and  approximate  temporal  averaging  are  equivalent  when  the  particle 
interarrival  time  becomes  long;  this  also  indicates  that  approximate 
temporal  averaging  does  not  approximate  true  temporal  averaging  when  the 
interarrival  time  is  long.  The  different  symbols  in  Fig.  20b  represent 
simulations  for  different  u,  and  /V  values  to  illustrate  that  this 
conclusion  is  not  limited  to  just  one  velocity-time  history. 

To  examine  the  arithmetic  and  approximate  time-integrated  mean 
velocity  for  turbulent  flow,  the  simulation  was  repeated  using  turbulent 
jet  hot-wire  measurements  as  the  velocity-time  function.  The  results, 
presented  in  Fig.  21 ,  again  show  no  bias  reduction  for  arithmetic  averages 
of  theoretical  IRLV  data  and  show  a  biasing  of  the  approximate  time- 
integrated  means  at  large  interarrival  times  corresponding  to  low  L^/At. 

As  in  the  sinusoidal  simulation,  a  significant  bias  reduction  of  arithmetic 
averages  could  only  be  effected  with  impractically  low  data  rates,  i.e., 
less  than  two  samples  per  hour. 

The  predicted  biasing  of  approximate  time-integrated  averages  of 
IRLV  data  is  particularly  important  if  approximate  time-integration 
is  used  to  avoid  statistical  biasing,  as  suggested  in  Refs.  20  and  48. 

As  Fig.  21  illustrates,  approximate  time-integration  will  produce 
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Figure  21.  Simulation  program  predictions  of  the  velocity  bias  of  a  turbulent  flow. 
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b.  Approximate  time-integrated  average,  v 
Figure  21.  Concluded. 
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averages  that  can  be  guaranteed  unbiased  only  if  the  data  rates  are 
high  enough  that  several  samples  are  recorded  per  integral  time  scale. 
This  would  require  data  rates  of  4,000  to  10,000  samples/sec  in  the 
experiments  conducted  in  this  investigation.  Most  IRLV  data  systems 
are  not  able  to  process  and  record  data  at  such  rates.  Also,  it  is 
often  not  feasible  to  seed  flows  heavily  enough  to  attain  the  necessary 
data  rates  without  possibly  altering  the  fluid  velocity  field. 

2.4.2  Data  Processing  Rate 

Statistical  bias  should  result  when  the  time  between  successive 

A 

IRLV  samples,  At.,  is  inversely  proportional  to  V^,  the  fluid  velocity 

average  evaluated  over  At^.  However,  the  time  between  successively 

recorded  IRLV  realizations,  At.,  may  not  be  inversely  proportional  to 
A  J 

V.  even  though  the  time  separating  each  particle  entrance  into  the 
probe  volume  obeys  Eq.  (9).  This  can  happen  simply  if  the  time  between 
some  particles  is  less  than  the  time  needed  to  evaluate  and  record  LV 
realizations.  For  example,  if  the  IRLV  instrumentation  can  only  record 
10  samples  per  second  but  an  average  of  1,000  particles  enters  the  probe 
volume  each  second,  then  a  particle  will  probably  be  in  the  volume  when 
the  LV  is  ready  to  process  another  sample.  Hence  the  time,  At_. ,  between 
recorded  realizations  will  be  nearly  constant  at  0.1  second  and  from  the 
discussion  following  Eq.  (2)  arithmetic  averages  will  be  unbiased. 


To  analyze  the  effect  that  processing  time,  k^,  can  have  on 
statistical  bias,  the  time  between  recorded  samples  is  expressed  as 


At.  -  k,  +  At' 

i  1  1 


(40) 


where  At*  is  the  delay  time  after  the  completion  of  a  data  acquisition 

1 

cycle  until  the  next  particle  arrival.  From  Eq.  (2)  the  mean  velocity 
can  be  expressed  as 

M 

(41) 

1=1 


E  V  Ck  I-  A  I*) 

j  t  j 

i=t  _ _ 
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<  V  >  1  + 
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If  the  processing  time  is  short  then  <At*>/k^  is  much  greater  than  unity 

and  At*  will  approach  At^  so  that  V  -  <V>  as  given  by  Eq,  (6) 

for  the  mean  velocity.  The  bias  will,  accordingly,  approach  that  previously 

* 

predicted.  However,  if  At  /k^  is  much  less  than  one,  as  when  the  processing 
time  is  long  compared  to  the  particle  interarrival  times,  the  approximate 
integral  mean  velocity  will  approach  the  arithmetic  mean  and  no  bias  should 
occur. 


To  estimate  the  relationship  between  the  temporal  and  arithmetic 
averages  of  the  data  for  intermediate  values  of  processing  time,  it  is 
necessary  to  estimate  the  form  of  the  correlation  coefficient 
From  Eq.  {6) , 

V  =  <  V  >  |l  +  RvAt| 

where,  by  evaluating  Eq.  (7)  for  the  mean  velocity, 

RvAt  =  <  V  '  At  '>/<  V  >  <  At  > 


(42) 

(43) 


Equation  (20)  requires  further  that 


V  =  <  V  > 


-i 


(44) 


Equating  Eq.  (42)  and  Eq.  (44)  shows  explicitly  that  the  velocity 

and  time  interval  fluctuations  are  anticorrelated  since 

2 


R 


VAl 


1  + 


(45) 


If  the  processing  time,  kfc,  is  short  compared  to  the  interval  between 
particle  arrivals  then,  from  Eq.  (40),  At^  -  At^  and  RyAt*  -  RyAt  This 
approximation  will  be  made  for  all  values  of  <At*>/kt<  Substituting  Eq. 
(45)  in  (41)  leads  directly  to 


V  =  <  V  > 


1  +  a—  +•  <Af> 

V  kt 

(1  +  1  +  <4ii> 

k  ' 


(46) 
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The  bias  parameter  corresponding  to  this  expression  is  shown  in  Fig. 

22  as  a  function  the  dimensionless  processing  delay,  kt/<At*>.  As  the 
delay  becomes  large,  the  arithmetic  and  temporal  averages  converge; 
the  net  effect  of  a  relatively  slow  data  acquisition  system  is,  then, 
to  reduce  bias.  Processing  delay  times  are  not  absolutely  constant  as 
assumed  in  the  preceding  analysis.  For  example,  the  delay  may  be  re¬ 
garded  as  the  sum  of  processor  recycle  time,  signal  processing  time, 
and  data  recording  time.  The  processor  recycle  time  is  approximately 
0.04  psec  for  the  Model  8  DDP  used  in  this  study.  The  data  acquisition 
system  can  record  data  at  rates  over  13,000  samples  per  second,  so  the 
time  associated  with  recording  was  at  most  77  psec.  The  signal  fre¬ 
quency  range  in  the  experiment  was  12.5  to  17  MHz  so  that,  when  the 
pulse  stretching  and  required  number  of  signal  cycles  associated  with 
the  DDP  were  considered,  signal  processing  times  ranged  from  60  to  80 
psec,  Consequently,  the  total  delay  between  successive  samples  was, 
at  most,  160  psec.  The  maximum  average  data  rates  encountered  were 
approximately  4,000  samples/sec,  making  kt/<At*>  *  0.7  in  the  worst 
case.  Figure  22  shows,  therefore,  that  the  bias  magnitudes  associated 
with  the  present  system  should  always  be  greater  than  60  percent  of 
the  maximum  possible  in  the  experiments.  Even  when  multidimensional 
effects  and  instrumentation  uncertainties  were  also  considered,  the 
present  experiments  should  have  detected  any  consistent  biasing  of  the 
IRLV  results,  but  none  was  found. 

The  random- samp ling  algorithm  (Ref.  25)  should  affect  bias  in 
much  the  same  way  as  a  constant  processing  delay,  except  that  the 
algorithm  introduces  random  rather  than  constant  delays  between  velocity 
realizations.  Thus,  since  constant  delays  are  effective  only  when  the 
delay  time  is  long,  random  sampling  will  be  effective  only  if  the  random 
delays  are,  on  the  average,  long  and  only  a  few  samples  are  selected  from 
a  large  data  set,  as  shown  by  Durao  and  Whitelaw  (Ref0  25).  To  verify 
the  effectiveness  of  the  procedures  used  to  randomize  the  experimental 
data,  the  same  procedures  were  used  to  randomize  data  generated  by  the 
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Figure  22.  Bias  reduction  resulting  from  minimum  processing  time. 
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numerical  simulation  with  the  hot-wire  velocity-time  function.  Random 
selection  of  500  samples  from  sets  of  2,000  often  effected  a  significant 
bias  reduction,  regardless  of  the  selection  rate,  but  only  when  the 
data  rate  in  the  numerical  simulation  was  high.  For  low  simulated  data 
rates  random  selection  was  often  ineffective.  Thus,  digital  randomization 
is  not  as  effective  as  the  random  time  processor  on-off  circuit  described 
in  Ref.  25  because  digital  randomization  cannot  be  expected  to  be 
equivalent  to  a  random  time-sampling  procedure.  However,  since  the 
digital  randomization  was  occasionally  effective,  it  should  provide 
evidence  of  significant  bias  of  IRLV  measurements  not  made  using  a  random 
time  sampling  circuit. 

2.4.3  Particle  Concentration  Fluctuations 

The  correlation  between  velocity  and  interarrival  time  fluctuations, 
which  forms  the  basis  for  statistical  bias  theories,  assumes  that 
particle  concentration  at  a  point  is  invariant.  This  condition  is 
extremely  hard  to  obtain  in  practice,  as  was  noted  in  Ref.  39  and 
observed  in  the  present  experiments.  As  an  aid  in  understanding  the 
experimental  results,  a  theory  can  be  developed  for  those  cases  where 
the  particle  concentration  as  well  as  fluid  velocity  is  time  dependent. 

For  this  purpose  a  formal  averaging  procedure  first  reported  by 
Cline  and  Bentley  (Ref.  41)  can  be  used.  The  fluid  velocity  at  a  point 
is  accordingly,  defined  by 

M 

*=i  (47) 

where  PCf^)  i-s  the  probability  that  the  variable  f  will  have  the 
particular  value  f  .  Likewise,  the  arithmetic  averages  of  the  individual 

K. 

realizations  of  velocity  are 


M 
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where  is  the  number  of  realizations  that  occur  in  some  interval 

f  +  Af  and  S  is  the  total  number  of  samples.  The  ratio  S,  /S  is  the 

k  M 

probability  that  the  particular  value  occurs  in  the  IRLV  data.  Now 

Sk  is  the  product  of  the  total  sampling  time,  T  ,  and  \  the  sampling 

rate  for  the  occurrence  of  the  f,  events.  Thus, 

k  ’ 

Sk  =  Vk  =  TtNkAVkp<V  (49) 

and  Eq ,  (48)  becomes 

M 

YWv 

<  f  >  =  _  (50) 

M 

^NkVkP<V 

k=  1 


The  particle  concentrations  and  velocities  corresponding  to  the 
k  interval  are,  respectively. 


Nk  =  Nk  +  K£ 


Vk  =  vk  +■  V' 


(51) 


so  that  Eq.  (50)  may  be  rewritten  for  the  mean  velocity  as 


<  V  >  =  V|  1  +  CTy  +  2Rnv  +  rnv21 


1  +  R 


(52) 


NV 


in  which  the  correlations  between  the  concentration  and  velocity 
fluctuations  are 


R  „  -  N'(V')n 

*1  M  V  11  —  - 


'NV  ' 


N  V 


(53) 


ru  2 


Similarly,  if  f  -  (V  -  V)  ,  Eq.  (50)  may  be  evaluated  for  the  turbulence 
intensity  in  a  flow  characterized  by  particle  concentration  fluctuations: 


'<v> 


1  +«^+2Rnv 


-  *7  -  RNV  +  So-  (1  -  Rnv)  -  Sv 

a— 


(54) 
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In  Eq.  (54)  third  order  or  higher  number-density  velocity  correlations 
have  been  neglected.  Note  further  that  jjhen  all  correlations,  n  =  0 
Eqs.  (52)  and  (54)  reduce  to  Eqs.  (20)  and  (22),  respectively,  thus 
establishing  the  equivalence  of  the  two  averaging  techniques. 

The  predicted  mean  velocity  and  turbulence  intensity  bias  for 
various  number  density-velocity  correlations  are  shown  in  Figs.  23  and 
24.  As  shown,  even  a  small  number  density-velocity  correlation  can 
significantly  alter  the  bias.  No  definitive  conclusions  can  be  reached 
as  to  the  sign  of  the  correlations  in  the  present  experiment. 

However,  the  R^v  correlation  can  be  expected  to  be  similar  to  the  species- 
velocity  correlation  in  flows  in  which  several  gases  are  mixing;  these 
correlations  are  generally  positive  (Ref.  49)  and  should  increase  the 
bias. 


Other  mechanisms,  such  as  particle  lag  effects,  may  reduce  the  bias 
in  IRLV  experiments.  However,  on  the  basis  of  the  preceding  analyses, 
where  none  of  the  postulated  bias-reducing  mechanisms  alone  was  adequate 
to  eliminate  bias  from  the  present  data,  it  is  unlikely  that  another 
mechanism  will  eliminate  the  bias.  A  possible  explanation  for  the  lack 
of  bias  in  the  present  data  may  be  simply  that  the  cumulative  effect  of 
several  mechanisms  reduced  the  bias  below  measurable  levels. 

3.0  SIGNAL  AMPLITUDE  STUDY 

Laser  velocimeters  can  provide  detailed  information  on  the  char¬ 
acteristics  of  a  flow  field  from  measurements  of  the  velocity  of  light¬ 
scattering  particles  contained  in  a  fluid.  In  most  applications, 
however,  the  fluid  velocity  and  not  the  particle  velocity  is  of  prime 
importance.  Therefore,  measurements  must  be  taken  to  ensure  that  particle 
response  to  fluid  transients  meets  a  desired  accuracy  criterion.  A 
satisfactory  approach,  where  feasible,  has  been  to  perform  particle  dynamic 
calculations  to  estimate  the  largest  particle  size  that  adequately  re¬ 
sponds  to  the  flow  and  then  seed  the  flow  with  particles  of  that  size  or 
smaller. 
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Figure  23.  Magnitude  of  mean  velocity  bias  for 
particle  concentration  fluctuations. 


Figure  24.  Magnitude  of  turbulent  intensity  bias  for 
fluctuations  in  particle  concentration. 
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Several  difficulties  are  associated  with  the  seeding  approach.  For 
example,  particle  response  caluclations  require,  a  pK-LafvL,  an  estimate 
of  the  fluid  velocity  and,  in  many  cases,  knowledge  of  the  turbulent 
structure  of  the  flow.  Lacking  this  basic  informtion,  one  can  only 
try  to  obtain  data  from  the  smallest  particles  detectable  by  the  LV 
used.  Seeding  the  flow  with  small  particles,  moreover,  does  not  guarantee 
that  all  of  the  data  obtained  arise  from  these  particles.  Agglomeration 
(see  Section  2.3)  or  the  presence  of  large  detectable  particles  in  the 
air  supply  may  introduce  measurements  that  make  up  a  significant  portion 
of  the  data  set.  Furthermore,  seeding  is  not  practical  or  desirable 
in  many  applications  (Refs.  5  and  50);  reliance  then  must  be  on  intrinsic 
particulates  whose  sizes  may  range  over  one  or  two  orders  of  magnitude. 

Many  of  these  difficulties  could  be  avoided  if  an  LV  could  be 
constrained  to  accept  only  data  from  particles  whose  sizes  are  less 
than  some  predetermined  value.  Prompted  by  these  considerations, 

Pfeiffer  and  von  Stein  (Ref,  36)  suggested  that  the  amplitude  of  a 
Doppler  burst  could  be  used  to  discriminate  against  large  particles. 

Data  obtained  downstream  of  an  oblique  shock  substantiated  this  claim. 

On  the  other  hand,  Andrews  and  Seifert  (Ref.  35)  had  only  modest  suc¬ 
cess  in  sizing  particles  on  the  basis  of  LV  signal  amplitudes.  Holve 
and  Self  (Ref.  51),  moreover,  observed  that  the  light-scattering 
amplitude  characteristics  of  particles  traversing  a  nonunif ormly 
illuminated  region  (such  as  the  probe  volume  of  an  LV)  are  trajectory 
dependent  so  that  particles  of  widely  varying  sizes  may  generate  bursts 
of  similar  amplitudes.  In  view  of  these  considerations,  the  concept  of 
discriminating  against  particles  on  the  basis  of  Doppler  burst  signal 
amplitude  presents  interpretation  problems  that  require  further  study 
before  adoption  into  LV  systems. 

3.1  MIE  SCATTERING  CHARACTERISTICS 

When  an  electromagnetic  wave  traverses  a  local  inhomogeneity  in 
refractive  index,  a  portion  of  the  wave  energy  is  redirected  (scattered) 
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as  the  result  of  reflection,  refraction,  and  diffraction  occurring  at 
the  wave-inhomogeneity  interface.  The  analysis  of  such  scattering 
phenomena  can  be  considerably  simplified  if  the  dimension  of  the  local 
disturbance  is  considerably  less  than  the  wavelength  of  the  electromagnetic 
radiation  (Rayleigh  scattering)  or,  conversely,  if  it  is  much  greater 
(geometric  optics  theory)  .  Unfortunately  disturbances  comparable  in  size 
to  particles  of  practical  interest  in  LV  measurements  of  fluid  ve¬ 
locity  are  generally  of  the  same  order  of  magnitude  as  the  wavelength. 
Therefore,  the  full  Mie  theory  must  be  considered  to  estimate  accurately 
the  magnitude  of  the  scattered  radiation.  Details  of  the  solution  of 
the  boundary  value  problem  solved  by  Mie  (Ref.  52)  and  simplifications 
to  the  theory  for  particles  have  diameters  much  different  from  the 
illuminating  wavelength  are  discussed  at  length  in  Refs.  53,  54,  and 
55.  The  theory,  consequently,  will  only  be  summari2ed  here  to  give 
sufficient  background  for  understanding  the  various  results  reported 
below. 

Consider  a  plane  electromagnetic  wave  travelling  in  an  infinite, 
homogeneous,  isotropic,  dielectric  medium.  If  the  electric  field,  2, 
and  magnetic  field,  H,  vectors  are  oriented  in  the  (x,y)  plane  of  the 

Cartesian  coordinate  system  shown  in  Fig.  25,  then  the  direction  of 

t  — y  ~y  — ^ 

propagation  given  by  the  Poynting  flux  vector,  S  =  E  x  H  is  directed 

along  the  z  axis.  The  intensity  of  light  at  a  point  P  =  P(r,S,4i) 

attributable  to  scattering  from  a  small  homogeneous  isotropic  sphere 

located  at  the  origin  of  the  coordinate  system  is 

J  =  ^_A_^  2  1  ■  F  {$  ,  4>)  (55) 

where  J  has  the  dimensions  of  (energy  time  1  area  -1 )  y  X  is  the  wavelength 
of  the  illuminating  wave,  I  is  the  intensity  of  the  illuminating  wave, 
and  F(0,$)  is  an  amplitude  function  describing  the  angular  distribution 
of  scattering.  For  computational  convenience  the  electric  field  vector 
can  be  resolved  into  components  perpendicular  and  parallel  to  the 
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"polarization"  plane  which  is  formed  by  the  incident  wave  and  the  radius 
directed  through  the  observation  point.  Accordingly,  the  intensity 
attributable  to  the  component  of  the  electric  field  vector  perpendicular 
to  the  scattering  plane  is  (Ref.  55) 

2 

I  S  .  I  sinz  th 

(56) 


(2ir) 


1  SJ  sin2  <6 


whereas  that  parallel  to  the  plane  is 

2 


J„ 


(£) 


I 


Sjj  |  cos2  (f) 


(57) 


The  amplitude  functions 
are  given  by  the  series 


required  for  evaluating  Eqs.  (56)  and  (57) 


Figure  25.  Mie  scattering  coordinate  system. 
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and 


S„  =  }  2n  + 1  i  A  r  ( cos  0)  -  B  w  (cos  0)1 
B  ^  n  (ti-1)  \  n  11  11  n  j 

n=l 

The  functions  t  and  tt  are  only  dependent  on  the  scattering 
n  n 

0,  and  are  defined  in  terms  of  the  Legendre  polynomials  of  a 
x*  by  the  relations 


= 


dP  (x*) 

n 

d  \* 


Tn 


x*  nn  (x*) 


-  (1  -  x-)2 

d  x* 


(59) 


angle, 

function 


(60) 

(61) 


In  evaluating  the  coefficients  of  the  amplitude  functions,  however,  the 
size  parameter,  a  =  7rd  /A,  the  ratio  of  the  particle  refractive  index 
relative  to  that  of  the  medium,  m  =  m  /m2,  and  the  product  6  =  ma  must 
be  introduced.  The  coefficients  may  then  be  expressed  as 


;i+T  |~Qnte)  QnQ3)  -  mQ>)  Qn(/3)~ 
Ll>n(aiQ'^)  -  0n(j8)_ 

^  |~Q;  fa)  Qn(fi)  -  mQn(q)  Q'JP)  ' 
|_  Qn  </3)  -  Q'l/S)  _ 


(62) 


(63) 


where  primes  denote  differentiation  of  a  function  with  respect  to  its 
argument.  The  Riccati-Bessel  functions,  Q^Cx*)  and  <t^(x*),  are  defined 
in  terms  of  the  Bessel  functions  of  the  first  kind  of  half-odd-integral 
order  through 


Q  „(**) 


Ik*) 


(64) 


and 


1>n(x*)  -  Qn(s+) 


(x+) 


(65) 
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With  the  above  relations,  local  values  can  be  computed  for  the 

intensity  and  phase  of  light  scattered  by  spherical  particles.  For 

large  particles,  the  series  given  by  Eqs.  (58)  and  (59)  converge 

slowly  so  that  many  terms  must  be  computed  to  minimize  truncation 

errors.  The  coefficients  A  and  B  ,  moreover,  are  sensitive  to  pre- 

n  n 

cision  errors  in  calculation  so  that  evaluation  of  the  scattering 
intensity  for  a  wide  range  of  parameters  is  practical  only  by  computer. 
In  recent  years,  several  programs  having  the  versatility  needed  for 
LV  studies  have  been  documented  (Refs.  56,  57,  and  58).  These  programs 
employ  recurrence  relations  for  the  Bessel  functions,  double  precision 
accuracy,  and  can  evaluate  the  light-scattering  amplitude  for  B  £  1,000. 
The  program  by  Dave  (Ref.  58),  as  modified  in  Ref.  59  to  account  for 
media  with  complex  refractive  indices,  was  used  for  the  Mie  scattering 
calculations  for  the  program  MIELV  described  in  Appendix  B. 

Several  definitions  are  useful  in  presenting  the  Mie  scattering 
calculations.  The  intensity  functions  are  defined  as 

ij_  1  (?)  =  S±  (0 i  -  (66) 

>„  <8)  =  |  S|(  (0)  |2  (67) 

For  either  6=0  deg  or  180  deg,  i  (6)  =  i  (6)  =  i(6)  so  that  there 

±  .1 

is  no  distinction  between  perpendicular  or  parallel  polarization  of  the 
light  wave  at  these  scattering  angles.  It  is  convenient,  therefore,  to 
nondinensionalize  the  intensity  function  with  i(0).  The  resulting 
normalized  intensity  functions  are  labeled  i*(6). 

The  influence  of  particle  size  on  the  on-axis  forward  (9=0)  and 
backscatter  (6  =  180  deg)  intensity  functions  is  shown  in  Fig.  26  for 
particles  scattering  light  in  a  nonabsorbing  medium  for  a  relative  re¬ 
fractive  index  of  a  =  1.33.  For  wavelengths  on  the  order  of  0.5  pm, 
such  as  those  used  in  the  LV  described  in  Section  2.2,  the  particle  size 
range  corresponding  to  Fig.  26  is  0.3  <  d^  <  16  pm.  The  relative 
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Figure  26.  Variation  of  Mie  intensity  with  particle  size  parameter 
for  water  droplets  (m  =  1.33). 
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refractive  index  corresponds  to  water  droplets  in  air.  The  forward- 
scatter  intensity  increases  nearly  monotonically  with  the  size  parameter, 
but  the  functions  for  backscatter  vary  in  an  irregular  manner  for  a  >  15. 
Generally,  the  intensity  functions  are  extremely  variable  in  any  of 
the  primary  variables  (6,  m,  or  a),  and  interpolation  between  computed 
values  is  subject  to  error. 


To  assess  the  influence  of  refractive  index  on  the  results,  the 
forward-scatter  intensity  function  for  several  particle  sizes  and 
refractive  indices  up  to  m  =  2.0  is  shown  in  Fig.  27.  Refractive  index 
clearly  has  a  pronounced  effect  on  the  intensity  function  for  small 
particles  where  i(9)  varies  over  two  orders  of  magnitude,  but  the  in¬ 
fluence  diminishes  with  increasing  particle  size.  Small  changes  in 
refractive  index,  however,  may  result  in  greater  variations  in  the 
intensity  function  for  large  particles  because  of  the  increased 
variability  in  i(0).  Representation  of  the  intensity  function  in 
backscatter  is  more  complex,  as  would  be  expected  from  Fig.  26  where 
it  is  shown  that  small  changes  in  a  lead  to  wide  variations  in  i(180 
deg).  For  convenience,  the  normalized  intensity,  i*(180  deg),  is 
plotted  versus  the  parameter  8  =  ma  to  minimize  overlapping  of  the 
results  with  increasing  refractive  index.  Typical  results  in  Fig.  28 
show  that  small  changes  in  refractive  index  as  well  as  particle  size 
can  lead  to  significant  variations  in  the  backscatter  intensity  function. 


The  effect  of  scattering  angle  on  the  normalized  intensity 

it  fi 

functions,  i^(0)  and  1^(6),  for  various  relative  refractive  indices 
is  shown  in  Figs.  29  and  30  for  particle  size  parameters  of  a  =  3.053 
and  12.212,  respectively.  The  value  m  =  1.46  is  typical  of  oleic  acid 
(Ref.  60),  and  m  =  1.72  is  representative  of  the  Al20^  particles  (Ref. 

61)  used  as  seed  material  in  the  bias  study.  For  7  =  0.5  pm,  the  particle 
sizes  represented  by  the  figures  are  nominally  d  “  0.5  and  2.0  pm.  For 
the  smaller  particle  (Fig,  29),  significant  decreases  in  scattering 
intensity  are  noted  near  60  and  150  deg  for  the  perpendicularly  polarized 
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Figure  27.  Effect  of  refractive  index  on  Mie  intensity 
function  in  forward  scatter  (0  =  0  deg). 
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results,  but  these  effects  are  less  pronounced  for  parallel  polarization. 

The  results  for  the  larger  particle  (Fig.  30)  show  that  the  situation  is 

decidedly  more  complex  at  large  values  of  a.  This  suggests  that  an  LV 

with  the  polarization  component  perpendicular  to  the  scattering  plane 

could  have  greatly  decreased  signal  amplitudes  for  some  particle  sizes 

and  refractive  indices.  Note  that  the  normalized  intensity  function 

for  on-axis  backscatter  detection  is  considerably  greater  (i*  -  10  '  to 
-2  -3 

10  )  than  the  i*  =  10  value  frequently  quoted  in  general  discussions 

(e.g..  Ref.  62)  of  LV  signal  amplitudes. 


3.1.1  Collector  Optics 


The  preceding  results  show  that  the  signal  amplitudes  to  be  expected 
from  particles  traversing  the  probe  volume  will  vary  significantly  with 
changes  in  particle  size,  refractive  index,  or  scattering  angle.  In 
practice,  however,  the  collector  optics  of  an  LV  will  integrate  the 
scattered  signals  over  a  region  in  space  which  will  tend  to  damp  the  ex- 
teme  variations  noted,  especially  in  backscatter.  The  rate  at  which  energy 
is  received  by  the  collector  lens  is 


EC  ■  #  JdSc 

s 

c 


(68) 


where  is  the  collector  lens  surface  area.  Working  with  the  polarization 
components  described  by  Eqs.  (56)  and  (57)  suggests  introducing  the  de¬ 
finitions 


E(.  =  *  $  r  2  ^1  sin2  dSc 


(69) 


E  =  /  jL\  I  $  r  2  i  (0)  cos2  d>  dS 

*  W  S  (70) 

C 

Although  the  integrations  indicated  In  Eqs.  (68)  through  (70)  are 
amenable  to  closed-form  solution  for  certain  collector  lens  geometries, 
more  general  results  can  be  obtained  by  numerical  techniques.  In  MIELV 
(Appendix  B) ,  the  integration  is  performed  for  a  lens  of  radius  R  of 
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focal  length  L  located  at  L  =  r  in  the  scattering  coordinate  system  of  Fig. 

2  2 

25.  The  solid  angle  of  the  collector  is  then  s2  =  ttR  /L  and  the  lens 
f-number  is  f  =  L/2R.  In  Figs.  31  and  32,  the  variation  of  the  dimen¬ 
sionless  captured  energy  with  scattering  angle  is  shown  for  a  range  of 
f-numbers.  For  these  calculations  both  circular  (i.e.,  random)  polarization 


and  a  refractive  index  of  m  =  1.54  corresponding  to  natural  particulates 
in  air  (Ref.  63)  have  been  used.  The  indicated  scattering  angle  corresponds 
to  the  midpoint  of  the  lens.  For  the  small-particle  (a  =  3.05)  and  large- 
particle  (a  =  12.2)  calculations  a  decrease  in  lens  f-number,  and  a 
corresponding  increase  in  the  lens  solid  angle,  smooths  local 

variations  in  the  detected  energy.  Consequently,  the  scattering  de¬ 
tected  by  the  LV  collector  optics  can  be  made  less  erratic  by  simply 
increasing  the  collector  lens  solid  angle.  This  trend  will  obviously 
enhance  the  ability  of  an  LV  system  to  discriminate  against  large 
particles  in  backscatter  systems  because  severe  decreases  in  signal 
amplitude  could  cause  the  signal  amplitude  of  a  large  particle  to  ap¬ 
proach  or  be  less  than  the  signal  level  of  a  smaller  particle  at  some 
scattering  angles. 

3.1.2  Illuminating  Intensity 

Because  of  the  nonuniform  illuminating  intensity  in  the  LV  probe 
volume,  signal  amplitudes  perceived  by  the  Doppler  data  processor  depend 
upon  the  trajectory  of  a  particle  as  it  sweeps  through  the  fringe  system. 

To  examine  the  trajectory  dependence  both  the  geometry  and  apparent 
variation  of  intensity  in  the  probe  volume  must  be  considered.  Both 
matters  have  been  discussed  in  detail  in  Refs.  29,  62,  and  64,  and 
the  model  proposed  by  Brayton  (Ref.  64)  will  be  used  in  the  following 
development . 
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Figure  32.  Effect  of  lens  f-number  on  scattering  functions  (a  =  12.21). 
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The  probe  volume  of  a  dual-scatter  LV  system  Is  formed  by  the 
intersection  of  two  laser  beams  which,  ideally,  are  of  equal  intensity, 
Gaussian,  and  similarly  polarized.  A  Cartesian  reference  frame  (X^, 

Y^,  Z^)  whose  origin  coincides  with  the  center  of  the  crossover  region 
is  depicted  in  Fig.  33.  The  bisector  of  the  illuminating  beams,  which 
defines  the  light  propagation  axis  of  the  LV,  is  directed  along  Z^,  the 
axis  is  coincident  with  the  detected  velocity  component,  and  X^  is 
normal  to  the  Y^Z  plane.  If  the  angle  between  the  beams,  0,  is  small, 
the  spatial  intensity  variation  is  given  by  (Ref.  64) 


(72) 


where  I  is  the  peak  intensity  in  each  of  the  Gaussian  beams  of  waist 
radius,  b,  which  form  the  probe  volume.  Since  only  a  fraction,  n»  of 
the  laser  power,  P  ,  actually  reaches  the  probe  volume,  the  peak 
intensity  for  a  one-component,  dual-scatter  LV  is 


*0  = 


»?pI 


ffb2  (73) 

Although  the  effective  probe  volume  is  limited  by  the  characteristics 

of  the  collector  optics  and  data  processor,  it  is  useful  to  define  a 

geometric  probe  volume  bounded  by  the  surface  where  the  signal  has 

2 

exponentially  decayed  to  1 /e  of  its  most  intense  value.  From  Eq. 

(72),  this  surface  is  described  by  the  ellipsoid 


(74) 


Thus  b  is  the  semiminor  axis  of  the  ellipsoid  and 


a  =  2k 

e 


(75) 


is  its  semimajor  axis.  A  set  of  nondimensional  coordinates,  (X*1  =  X  /b, 
Y^*  =  Y^/b,  Zn*  =  Z^/a)  may  be  introduced  which  transform  the  geometric 
probe  volume  to  a  sphere  and  allow  Eq .  (72)  to  be  rewritten  as 


[-2  (X’j 2  +  Y 


+  2 


•;2>]  j 


cosh  4  Y  J  Zj  +  LOS  7T  Nj  Y 


(76) 
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Figure  33.  Probe  volume  coordinate  system. 


where  I*  =  I/I  and  Nr  =  2b/K.  is  the  number  of  fringes  in  the  plane 

o  f  o 

Z.j*  =  0.  The  maximum  probe  volume  intensity  occurs  at  the  origin  of 
the  coordinate  system  where  I  =  21^.  The  fringe  model  [Eq.  (76) J 
of  the  probe  volume,  however,  requires  the  intensity  to  approach  zero 
in  an  alternating  fashion  as  depicted  in  Fig.  34,  where  the  intensity 
distribution  expected  for  the  LV  of  Section  2.2,1  is  shown  at  =  0 

for  three  values  of  Z^*.  The  maximum  Hie  scattering  amplitude  in  a 
dual-scatter  system  is  highly  dependent  upon  the  trajectory  of  the 
particle  as  it  traverses  the  probe  volume.  Even  though  this  dependence 
does  not  absolutely  preclude  discriminating  against  large  particles  on 
the  basis  of  their  signal  amplitudes,  such  discrimination  can  only  be 
accomplished  statistically  and  some  large-particle  data  will  generally 
be  present  in  any  data  set  obtained. 
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To  characterize  the  illuminating  intensity  without  treating  the 
details  of  the  interference  pattern  along  the  particle  trajectory,  several 
representative  illuminating  intensities  can  be  defined.  The  maximum 
intensity,  I  =  21^,  is  clearly  one  possibility.  A  second  is  the  average 
intensity  in  the  probe  volume 


] 


d  Z 


* 

i 


(77) 


where  I*  is  defined  by  Eq.  (76).  Note,  however,  that  a  processable 
signal  will  likely  result  from  particles  which  cross  the  central 
portion  of  the  probe  volume.  For  the  assumptions  that  (1)  the  flow 
is  one-dimensional  and  (2)  the  probe  volume  is  small  enough  that 
velocity  and  particle  concentration  are  uniform  across  it,  the  prob¬ 
ability  density  function  for  a  particle  trajectory  along  any  plane 
Z^*  =  constant  is 

f(Z*)  =  2  [l  -  Z*2 

1  77  \| 


(78) 


A  most  probable  Intensity  can  then  be  defined  by 

T  21  f 

71  J 

7  * 

L  l 


r*  •  kz*)  dz* 


(79) 


where  I*(Z*)  is  the  average  intensity  along  the  constant  Z*  planes, 
That  is. 


r  (z*)  = 


i 


n-  Cl  -  7. 


-Jf 

Y  *  V  * 
l  1 


I*  dX*  dY* 


(80) 


The  above  integrations  are  performed  numerically  by  the  subroutine 
PVIP  described  in  Appendix  C.  Typical  results  are  presented  in  Fig.  35, 
where  the  variation  of  I*  with  the  probe  volume  axial  coordinate  is  shown 
along  with  the  most  probable  and  average  intensities.  <I>  will  be  ap¬ 
proximated  as  21^/5.  Although  <I>  should  provide  representative  intensities 
for  a  dual-scatter  LV,  in  discriminating  against  large  particles  the 
maximum  signal  amplitude  is  of  primary  concern.  By  assuming  each  particle 
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will  cross  X*  =  Y*  =  0  while  traversing  the  probe  volume,  Eq.  (76)  reduces 

2 

to  I*  =  2  exp  (-2  Z*  )  and  a  most  probable  maximum  intensity  can  be 
max  1  , 

defined  as 


<Imax>  =  ^  j  -P  (-2Z*2)  -  Z) 2  dZ;  (81) 

O 

A  numerical  integration  of  Eq.  (81)  gives  <1  >  -  21  /3. 

to  max  D 

The  adequacy  of  these  representative  illuminating  intensities  will 
be  discussed  in  the  data  comparisons  later  in  this  section.  By  using 
a  single  representative  intensity,  the  Mie  scattering  analysis  is 
simplified  and  results  obtained  for  one  intensity  value,  I.,  are  readily 
scaled  to  another  value,  I^,  since,  from  Eqs.  (69)  or  (70),  I^/lj 

The  representative  intensities  make  it  possible,  moreover,  to  assess^- 
whether  the  LV  used  in  the  bias  study  can  be  unambiguously  adapted  for 
particle  size  discrimination.  Using  <Imax>»  the  effect  of  viewing  angle 
on  the  scattered  light  levels  is  shown  for  several  forward-scatter 
collection  angles  in  Fig.  36,  and  backscatter  results  are  given  in  Fig.  37. 
For  particle  diameters  up  to  10  um,  the  theoretical  signal  amplitudes 
increase  nearly  monotonically  at  several  angles  in  both  forward-  and  back¬ 
scatter  (particularly  0  =  40,  60,  and  170  deg).  Fortuitously,  a  0  =  170  deg 
collector  configuration  approaches  the  design  condition  for  the  LV  and 
traverse  system  of  Fig.  6. 
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Figure  36.  Effect  of  particle  size  on  light  scattered 
for  forward-scatter  system. 
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1  =  <^max> 

D  =  700  w/cm^ 
c  =  O.Olir  sr 
n  =  (1.54,  0) 
i  *  l/2(  1  j_  +  1„) 
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3.2  SIGNAL  AMPLIFICATION 

The  theoretical  signal  levels  in  Figs.  36  and  37  represent  the  rate 

at  which  scattered  light  energy  is  intercepted  by  the  collector  optics. 

To  discriminate  against  particles  whose  size  exceeds  a  predetermined 

value  requires  consideration  of  the  amplitude  of  the  electrical  signal 

at  the  Doppler  data  processor  (DDP) .  To  proceed  from  the  collected 

radiation,  E  ,  to  the  detected  voltage,  e  ,  consider  the  amplification 
c  a 

and  attenuation  of  the  signal  by  the  photomultiplier  tube  and  the  LV 
electronics.  For  this  study,  the  signal  level  will  be  related  to 
collected  energy  by  the  linear  relation 

(JU)(!.t0  <82) 

where  is  the  ratio  of  the  signal  level  after  amplification  to  its 

value  entering  the  i1"*1  amplifier  and  Tj  is  the  corresponding  ratio 

r  ,  th  J 

for  the  j  attenuation  process. 

For  the  specific  LV  being  studied,  only  gains  attributable  to 
photomultiplier  tube  (PMT)  amplification  need  be  considered  since  no 
other  amplifiers  were  used.  Similarly,  attenuation  was  primarily 
attributable  to  a  laser  line  filter  between  the  aperture  and  PMT,  as 

O 

shown  in  Fig.  1.  For  the  30-A  bandwidth  filter  used,  the  transmission 
ratio  at  A  =  0.5145  pm  was  =  0.765. 

A  photomultiplier  converts  incident  radiation  into  electrical  energy 
by  a  photoemission  process  and  then  amplifies  the  signal  by  means  of 
secondary  emission.  The  elements  of  a  typical  tube  may  be  regarded  as 
(1)  a  semitransparent  photocathode,  which  emits  electrons  because  of 
interaction  of  the  incident  radiant  energy  with  electrons  in  the  photo¬ 
cathode  material,  (2)  several  dynode  stages,  in  which  secondary  emission 
results  from  the  impact  of  electrons  on  the  dynode  surface,  and  (3)  an 
anode,  which  collects  the  electrons  emitted  in  the  last  dynode  stage 
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and  provides  the  output  current  for  the  PMT.  Although  photomultiplier 
tubes  can  be  analyzed  in  detail  by  quantum  statistical  techniques,  in 
the  present  application  a  PMT  will  be  regarded  simply  as  a  linear  device 
characterized  by  a  transfer  function  H  so  that 

Lj 


'A  = 

The  gain  attributable  to  the  PMT  is 


»G  r'< 


,  then 


V, 


l 


t: 

c 


T 


l 


(83) 


(84) 


where  R  is  the  LV  electronics  circuit  resistance.  The  gain  associated 
e  ° 

with  a  particular  photomultiplier  is,  however,  limited  by  its  supply 
voltage,  e  ,  and  dynode  voltage  divider  circuitry.  Typically,  H  for 

s  ^  c 

various  tubes  may  range  from  10  to  10  ma/yw  by  proper  tube  choice 
the  LV  signals  can  be  amplified  to  a  detectable  level,  A  serious  concern 
for  discrimination  studies  is  the  observed  variation  in  gain  for 
particular  tube  types  under  ostensibly  similar  conditions.  Variations 
in  amplification  levels  ranging  over  two  orders  of  magnitude  have  been 
observed  (Ref.  65).  Consequently,  discriminating  against  particles  on 
the  basis  of  signal  amplitude  must  be  based  on  calibration  of  the  PMT 
system  used  rather  than  on  manufacturer's  specifications. 


Such  calibrations  were  performed  for  several  photomultipliers 
using  the  arrangement  illustrated  in  Fig.  38.  In  this  calibration  a 
partially  transmitting  mirror  is  placed  at  the  center  of  the  probe 
volume  so  that  an  essentially  steady  anode  current  results.  Note  that 
the  laser  radiation  from  the  illuminating  optics  must  be  attenuated 
to  prevent  damage  to  the  PMT,  Then,  by  choosing  an  appropriate  re¬ 
sistance,  a  voltage  level  can  be  read  on  either  a  digital  voltmeter 
or  an  oscilloscope.  The  illuminating  intensity  was  varied  by  changing 
the  neutral  density  filters  used  to  attentuate  the  laser  radiation  and/ 
or  adjusting  the  laser  power.  After  each  series  of  measurements  the 
PMT  was  replaced  by  a  calibrated  laser  power  meter  so  that  the 
illuminating  intensity  could  be  measured. 
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Figure  38.  PMT  calibration  arrangement. 


Typical  calibration  curves  for  an  RCA  931A  PUT  coupled  with  a  voltage 
divider  circuit,  which  limited  the  anode  current  to  about  1  ma,  are  shown 
for  various  supply  voltages  in  Fig.  39.  If  the  amplification  function  is 


defined  as 


C 


(85) 


the  variation  of  H„  with  supply  voltage  can  be  determined.  A  comparison 

G 

of  the  experimental  values  with  the  manufacturer's  data  (Ref.  65)  is 
made  in  Fig.  40  and  emphasizes  the  importance  of  calibration  of  the 
particular  tube  in  adequately  predicting  signal  amplitude.  Also 
apparent  in  Fig.  39,  the  anode  current  varies  nonlinearly  with 
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Figure  39.  PMT  amplification  for  several  power  supply  voltages. 
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illuminating  intensity  once  PMT  saturation  is  approached.  Therefore, 
the  level  of  the  scattered  light  must  be  limited  if  the  linear  relation 
proposed  in  Eq.  (83)  is  to  hold.  The  observed  illuminating  intensity 
needed  to  produce  saturation  of  the  PMT  is  shown  in  Fig,  41  as  a  function 
of  supply  voltage.  Note  that  these  results  apply  to  a  particular 
PMT  and  cannot  be  generalized  to  include  all  PMT’s  of  this  type. 

3.3  EXPERIMENTAL  SIGNAL  LEVELS 

As  proposed  in  Section  3.1,  the  Mie  scattering  characterisitics 

of  particles  traversing  the  probe  volume  of  a  dual— scatter  LV  were 

estimated  by  considering  the  intersecting  beams  as  a  single  beam  at  some 

equivalent  illuminating  intensity.  The  most  probable  intensity, 

<I>  [Eq.  (79)]  ,  was  used  to  obtain  representative  signal  amplitudes,  and 

the  maximum  amplitude  of  the  signal  was  determined  from  the  probability 

weighted  maximum  intensity,  <1  >  (Eq.  81).  The  adequacy  of  these 

e  y  max 

equivalent  intensities  was  checked  by  experimentally  determining  LV 
signal  amplitudes  for  monodisperse,  spherical  particles  of  various  sizes. 
The  confined  jet  described  in  Section  2,2  provided  a  flow  free  of  in¬ 
trinsic  particles  for  nominal  exit  velocities  of  =  15  and  30.5  m/sec. 

A  Berglund-Liu  aerosol  generator  (Ref.  40)  provided  spherical,  oleic  acid 
particles  having  diameters  of  either  d^  =  2.0,  5.0,  or  10.0  am. 

Microscopic  sizing  of  collected  samples  verified  that  the  particles  were 
essentially  spherical  and  were  uniformily  of  the  nominal  diameters. 

Signals  were  generated  by  the  Bragg-diffracted  LV  operated  in  a  one- 
component  mode  at  a  collector  lens  position  of  0  =  168.4  deg.  At  this 
collector  angle,  stray  reflections  from  the  cell  and  optical  components 
were  negligible,  and  the  laser  line  filter  described  earlier  eliminated 
signals  from  radiation  sources  having  wavelengths  other  than  a  =  0.5145  pm. 
PMT  calibrations  showed  that  the  photomultiplier  transfer  function,  H  , 

i  G 

was  4  x  10  ma/pw  and  no  other  signal  amplification  was  used.  The  PMT 
anode  signal  was  terminated  with  an  =  50-ohm  impedence,  and  laser  power 
was  varied  with  particle  size  to  ensure  that  PMT  saturation  current  was 
not  approached. 
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To  determine  the  maximum  amplitude  of  signals  generated  by  particles 
traversing  the  probe  volume,  the  oscilloscope  trigger  was  calibrated 
using  a  pulse  generator  with  output  signals  of  known  amplitude  and  a  given 
repetition  rate.  This  signal  was  input  to  the  oscilloscope,  and  a  frequency 
counter  counted  the  oscilloscope  trigger  pulses  that  resulted  when  the 
calibration  signal  exceeded  the  trigger  level.  The  oscilloscope  trigger 
was  set  at  a  level  that  decreased  the  trigger  pulse  rate  to  half  the  pulse 
generator  output  rate  so  that  the  generator  output  amplitude  provided  a 
precise  setting  of  the  voltage  level  at  which  the  oscilloscope  would  trig¬ 
ger.  The  PMT  signal  was  then  input  to  the  oscilloscope  and  sample  rate 
measured  by  the  frequency  counter.  Data  were  taken  in  four  10-sec  intervals 
to  obtain  the  mean  data  rate  over  a  40-sec  interval.  The  process  was  then 
repeated  to  obtain  the  sample  rate  at  successively  higher  trigger  levels. 
Although  the  trigger  could  be  calibrated  in  2-mv  intervals,  10-mv  steps 
were  used. 

Fdt  the  PMT  and  voltage  divider  used,  anode  saturation  was  approached 
near  a  signal  level  of  50  mv.  Consequently,  the  laser  power  was  decreased 
as  the  particle  size  increased  to  avoid  the  nonlinear  PMT  amplification 
observed  (Fig.  39)  near  PMT  saturation.  Since  the  available  amplifiers 
caused  serious  signal  distortion  and  inversion,  the  unamplified,  low  PMT 
signal  levels  were  used.  At  these  low  signal  levels,  random  low-level 
noise  signals  were  detected  at  trigger  levels  below  9  mv.  This  was 
evidenced  by  observed  data  rates  that  far  exceeded  the  ability  of  the 
seeder  to  deliver  particles  to  the  probe  volume.  Therefore,  detection 
was  limited  to  signals  having  amplitudes  greater  than  10  mv,  and  the 
number  of  particle  signals  below  that  level  were  estimated  by  the  difference 
between  the  observed  and  expected  events. 

Despite  this  limitation,  remarkably  consistent  results  were  obtained 
in  all  tests.  Figure  42  shows  these  results  in  a  histogram  of  oscilloscope 
sampling  probabilities  corresponding  to  each  trigger  setting.  The  typical 
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signal  level  should  be  between  10.7  and  13.3  mv  in  each  case,  as  indicated 
using  the  computer  model  MIELV  with  an  illuminating  intensity  taken  as  the 
peak  single  line  intensity,  Iq,  for  each  laser  power  used.  Since  Fig.  42 
shows  that  the  most  probable  signal  levels  are  below  10  mv,  I  is  too  high  a 
value  to  be  used  for  the  illuminating  intensity  in  MIELV  calculations. 
Fortunately,  the  signal  level,  ed>  is  a  linear  function  of  the  collected 
energy,  E  ,  and  thus,  in  turn,  a  linear  function  of  the  illuminating 
intensity,  so  the  MIELV  predicted  signal  levels  will  agree  with  the  most 
probable  observed  signal  levels  if  the  illuminating  intensity  is  taken 
as  I  / 2  to  31  /4 .  This  compares  favorably  with  the  value  <1  >  =  0,67  I 

Q  Q  ITlclX  O 

predicted  by  Eq.  (81).  With  the  uncertainty  in  the  experiments,  however, 
this  result  can  only  be  taken  as  evidence  that  with  adequate  calibration 
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the  single-beam  model  of  an  LV  can  produce  signal  amplitudes  that  are 
representative  of  experimental  results.  Therefore,  program  MIELV  was  used 
to  assess  the  feasibility  of  sizing  or  discriminating  against  particles 
on  the  basis  of  LV  signal  levels. 

3.4  PARTICLE  SIZING  AND  DISCRIMINATION 

Figure  43  shows  the  variation  in  theoretical  signal  amplitude  for 

the  LV  used  in  these  studies  at  a  collector  angle  of  0  =  168.4  deg  and  an 

illuminating  intensity  corresponding  to  a  laser  power  of  P  =  150  raw. 

Jj 

The  particles  are  assumed  to  be  oleic  acid  in  air.  For  particle  diameters 

up  to  2.0  jm  the  signal  level  is  a  monotonically  increasing  function  of 

particle  size.  The  decrease  in  signal  level  for  particles  slightly 

greater  than  2.0  um,  however,  limits  the  useful  range  for  sizing  studies 

to  dp  _<  1.8  pm.  Therefore,  if  the  signals  were  suitably  amplified  the 

present  LV  could  size  oleic  acid  particles  in  the  range  0.5  <  d  >  1.8  m. 

P 

The  lower  limit  is  set  by  the  noise  associated  with  the  system.  For 

dp  2.0  pm,  however,  the  variations  in  signal  levels  with  changes  in 

particle  size  would  only  permit  gross  size  determinations.  That  is, 

while  a  2~um  particle  could  be  distinguished  from  one  whose  diameter  was 

5  pm,  it  would  not  be  possible  to  distinguish  between  3.0-  and  3.8-pm 

particles.  It  should  be  noted  that  it  is  the  particle  size  parameter, 

a  =  md  /A,  that  actually  determines  the  extent  of  the  region  where 

monotonic  increases  in  signal  amplitudes  are  expected.  Therefore,  if 

the  LV  system  used  a  laser  which  radiated  in  the  infrared  the  useful 

sizing  range  could  be  extended  as  noted  in  Ref.  35.  Specifically,  for 

a  wavelength.  A,  of  10  tint  the  monotonic  range  would  extend  to  d  <34  pm 

P 

if  all  other  LV  system  parameters  could  be  maintained. 

Although  the  use  of  signal  amplitudes  for  particle  sizing  is  limited 
to  specific  size  range,  the  use  of  signal  amplitudes  for  size  discrimination 
is  feasible.  For  the  conditions  applicable  to  Fig.  43,  discriminating 
against  all  signals  whose  amplitudes  exceed  =  8  mv  would  exclude  all 
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Figure  43.  Effect  of  particle  size  on  signal  amplitude. 
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particles  having  diameters  greater  than  2.0  pm.  The  variation  in  signal 
level  with  laser  power  or,  equivalently,  illuminating  intensity  is  shown 
in  Fig.  44.  For  a  laser  power  of  300  mw,  constraining  the  LV  to  process 
signals  of  amplitudes  less  than  e^  =  30  mw  would  eliminate  particles  with 
diameters  greater  than  3.25  pm  (as  well  as  some  particles  with  diameters 
slightly  less  than  3.25  pm).  The  figure  illustrates,  moreover,  the 
effect  of  setting  a  lower  threshold  level  either  deliberately  or  because 
of  processor  limitations  or  the  presence  of  low-level  noise.  For  the  above 
example  no  particle  below  =  1.75  pm  could  be  detected,  although  this 
situation  could  be  remedied  by  introducing  suitable  amplification  after 
the  PMT.  In  this  regard,  an  amplifier  with  a  gain  of  10  would  allow 
detection  of  1-gm  particles  and  cause  the  signals  below  30  mv  to  fall 
within  the  monotonic  range  of  Fig.  43,  thereby  eliminating  ambiguities  in 
discrimination  limits. 

The  concept  of  discriminating  against  large  particles  is,  perhaps, 
of  most  interest  in  a  flow  where  several  types  of  particles  are  present. 
Such  a  flow  would  commonly  occur  when  artifical  seeding  is  used  to 
enhance  data  rates  or  control  the  intrinsic  particle  characteristics.  The 
range  of  amplitudes  expected  for  air  containing  natural  particles 
(m  =  1.54),  Ref.  63,  aluminum  oxide  seed  particles  (m  =  1.72),  and  water 
droplets  (m  =  1.33)  is  shown  in  Fig,  45.  The  observations  made  on  the 
possibility  of  sizing  particles  apply  in  each  of  these  cases  as  well, 
although  the  region  where  the  amplitude  function  exhibits  monotonic 
behavior  is  a  strong  function  of  refractive  index.  Specifically,  for 
water  droplets  the  range  for  which  sizing  is  possible  is  <_  1 . 5  pm, 
whereas  for  either  A1?0^  or  intrinsic  scatters  the  limit  is  approximately 
dp  <_  0.8  pm.  In  Fig.  46,  these  data  are  replotted  for  various  illuminating 
intensities.  Again  the  LV  can  be  constrained  to  exclude  signals  from  large 
particles.  However,  low-level  discrimination  as  well  as  an  upper  limit 
must  be  used  if  some  large-particle  signals  (H^O,  d  “  3.0)  are  to  be 
eliminated  from  consideration. 
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Figure  44.  Signal  level  variation  with  illuminating  intensity. 


110 


Theoretical  Signal  Level)  e.,  mv 


A  E  DC  -T  R  -78-65 


1 _ I _ I _ I _ I _ 1 _ I _ 'III 

0  20  40  60  80  100 

2 

Illuminating  Intensity,  I0,  w/cm 


Figure  46.  Effect  of  illuminating  intensity  on  signal  amplitudes 
for  several  particulates. 
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In  summary,  the  use  of  the  Mle  scattering  amplitude  of  LV  signals  to 
size  particles  implies  a  degree  of  control  over  the  experiment  that  is  not 
always  practical.  Specifically,  use  of  an  LV  for  this  purpose  requires: 

1 .  the  freedom  to  choose  an  optimum  collector  angle  and  collector 
lens  so  that  a  monotonic  amplitude  function  can  be  obtained, 

l 

2.  knowledge  of  the  refractive  index  of  the  particles, 

3.  calibration  of  PMT  amplification  characteristics, 

4.  choice  of  laser  power  to  keep  the  PMT  from  saturation  yet 
provide  adequate  S/N  ratios  for  velocity .measurement,  and 

5.  low-noise  levels  in  electronic  signals. 

While  some  of  these  constraints  may  be  relaxed  in  discriminating 
against  large  particles,  some  knowledge  of  signal  amplitude  versus  particle 
size  must  be  obtained  by  calibration  or  calculation  to  set  the  discriminator 
voltage  properly.  Furthermore,  in  a  seeded  flow  with  intrinsic  particles 
a  low-level  discriminator  may  also  be  needed  to  eliminate  some  large- 
particle  signals.  In  addition,  the  results  are  extremely  sensitive  to 
each  of  the  parameters  characterizing  the  Mie  scattering  process.  As 
a  result,  extrapolating  results  to  condition  that  differ  even  slightly 
from  those  considered  theoretically  or  observed  experimentally  is 
unreliable.  Finally,  it  must  be  realized  that  particle  sizing  or 
discrimination  from  signal  amplitudes  can  only  be  accomplished  in  a 
statistical  sense  since  detailed  knowledge  of  where  each  particle  crosses 
the  probe  volume  is  unattainable. 
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4.0  CONCLUSIONS 

Analytical  and  experimental  studies  were  performed  to  determine 
(1)  whether  statistical  bias  in  the  mean  and  central  moments  of  ve¬ 
locity  determined  from  individual  realization  laser  velocimeter  (IRLV) 
measurements  is  a  significant  source  of  error  and  (2)  whether  the  maximum 
amplitude  of  LV  burst  signals  could  be  used  to  size  particles  or 
discriminate  against  signals  from  large  particles.  Both  problems  are 
important  in  the  interpretation  of  LV  measurements  since  affirmative 
answers  to  either  question  could  alter  both  the  course  of  Doppler  data 
processor  (DDP)  development  and  the  experimental  procedures  and  data 
reduction  techniques  used  in  LV-supported  test  programs. 

The  analytical  studies  of  statistical  bias  showed  that,  for  a 
bias  to  exist  in  IRLV  data,  the  fluid  velocity  fluctuations  and  the 
times  between  particle  arrivals  must  be  correlated.  While  the  sampling 
mechanisms  peculiar  to  IRLV  utilization  predict  such  correlations,  it 
was  shown  by  both  limiting  statistical  arguments  and  a  bias  simulation 
computer  program  that  the  correlation  function  is  altered  by  (1)  three- 
dimensional  velocity  fluctuations,  (2)  fluctuations  in  particle  concen¬ 
tration,  and  (3)  data  acquisition  rate  limitations  imposed  by  the  DDP  and 
data  acquisition  systems.  The  simulation  program  showed,  moreover,  that 
various  bias-correcting  algorithms  proposed  for  improving  IRLV  averages 
in  statistically  biased  situations  are  effective  only  at  high  data  rates 
(in  excess  of  two  samples  per  integral  time  scale)  and  might  fail  at  low 
data  rates. 

In  the  experimental  program,  mean  velocities  and  turbulence 
intensities  in  the  near  field  of  a  subsonic  jet  were  determined  (as 
applicable  from  pitot  pressures,  a  hot-wire  anemometer,  and  an  IRLV 
system) .  Although  care  was  taken  to  produce  conditions  favorable  to 
the  attainment  of  biased  results,  no  consistent  experimental  evidence 
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of  bias  was  found.  Furthermore,  random  selection  of  the  data  produced 
results  which  were  nearly  identical  to  those  obtained  from  the  original 
data  sets.  Consequently,  the  velocity-arrival  time  correlation  necessary 
for  obtaining  bias  was  negligible  and  no  bias  occurred. 

In  conclusion,  (1)  statistical  bias  is  not  present  in  all  IRLV 
data,  (2)  routine  correction  for  bias  can  produce  serious  errors,  and 
(3)  modification  of  DDP's  to  incorporate  bias-correcting  algorithms  or 
circuitry  is  not  always  warranted.  Since  random  sampling  can  be  readily 
accomplished  using  standard  computer  subroutines,  data  reduction  pro¬ 
cedures  should  include  a  check  for  bias  using  random  selection  of  the 
recorded  data  sets.  Bias-correcting  algorithms  should  be  used  only  if 
this  technique  shows  significant  bias.  Numerical  random  selection  should 
not  be  used,  however,  to  correct  for  bias. 

To  study  the  ability  of  signal  amplitude  to  discriminate  particle 
size,  a  computer  program  was  developed  to  (1)  determine  the  Mie  scattering 
amplitudes  of  signals  from  particles  traversing  the  probe  volume  of  a 
dual-scatter  LV  and  (2)  compute  the  signal  level  expected  at  the  data 
processor  after  amplification  and  attenuation  of  the  scattered  signal 
by  the  system.  For  simplicity,  the  program  (MIELV)  assumed  that  the 
illuminating  intensity  in  the  probe  volume  of  a  dual-scatter  LV  system 
was  equivalent  to  a  single-beam  intensity  determined  in  an  auxiliary 
program  (PVIP) . 

Experiments  were  performed  to  validate  the  equivalent  beam  assumption 
and  to  ascertain  how  closely  typical  photomultiplier  tubes  correspond  to 
manufacturer's  specifications.  It  was  determined  that  using  an  equivalent 
intensity  qualitatively  gives  signal  levels  that  conform  to  those 
experimentally  observed.  Photomultiplier  calibrations,  however,  showed 
an  order  of  magnitude  variation  in  the  amplification  obtained  with 
different  tubes.  For  this  reason,  any  sizing  or  discrimination  study 
must  be  performed  with  an  accurately  calibrated  tube. 
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Based  on  the  results  of  these  experiments,  (1)  the  wide  variation 
in  signal  levels  with  particle  size,  index  of  refraction,  and  collector 
optics  design  makes  it  imperative  that  sizing  or  discrimination  studies 
be  supported  with  detailed  Mie  scattering  calibrations,  (2)  the  predicted 
LV  signal  amplitudes  depend  upon  the  particle  trajectory  through  the  probe 
volume  so  that  particle  identification  can  only  be  performed  statistically 
or  by  restricting  the  collector  optics  to  view  only  a  portion  of  the 
probe  volume,  (3)  particle  sizing  is  possible  only  for  a  limited  range  of 
sizes  or  at  particular  collector  optics  scattering  angles,  (4)  the  ef¬ 
fective  sizing  range  increases  linearly  with  wavelength  and  can,  there¬ 
fore,  be  extended  by  utilizing  IR  light  sources,  and  (5)  discrimination 
against  large  particles  can  be  achieved  qualitatively  even  in  flows 
containing  several  types  of  particles. 

Future  DDP's  should  have  low-  and  high-level  discriminators  to 
identify  or  eliminate  signals  from  large  particles  during  signal 
processing.  The  computer  programs  developed  during  the  study,  moreover, 
should  be  valuable  in  the  design  of  LV  systems  for  applications  dif¬ 
fering  appreciably  from  those  considered  in  this  study.  Finally,  while 
beyond  the  scope  of  this  study,  the  single-beam  assumption  is  not 
essential  to  the  MIELV  program  and  should  be  removed  in  future  program 
developments . 
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APPENDIX  A 

RELATIONSHIP  BETWEEN  LONG-  AND  SHORT-TERM  AVERAGES 


-  In  pulsating  flows  the  temporal  mean  velocity,  V^,  over  time  interval, 

At  ,  may  differ  from  the  temporal  mean  velocity,  V,  defined  over  an 
1 

infinitely  long  time  interval.  The  difference  between  the  mean  veloc¬ 
ities  is  given  by 


S 

o 


The  variance  or  mean-square  value  of  this  difference  is 


K  +T 

n 


Q 


1  1ID 

T 


L 

rI 


V) 


(A-2) 


Assuming  that  the  velocity,  V(t),  is  statistically  stationary,  rearrang¬ 
ing  the  order  of  integration,  and  performing  a  Reynolds  decomposition  of 
the  velocity,  gives 


where  v(t)  is  the  velocity  fluctuation  about  V.  The  integration  in  the 
(t,  t')  plane  is  over  a  square  with  limits  0  and  t.  Since  the  integrand 
is  symmetrical  in  t  and  t' ,  the  double  integration  can  be  replaced  with 
twice  the  single  integration  over  the  triangle  that  is  half  the  square. 
This  leads  to, 
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(A-4) 


s 


To  alter  the  form  of  the  integrand  in  Eq.  (A-4)  the  difference, 
s,  is  defined  as  t,  so  that  Eq.  (A-4)  can  be  written 


The  integrand  is  now  recognized  as  the  velocity  autocorrelation,  defined 
in  coefficient  form  as 


RM 


l 


A 


]  i  m 

T 


T 

I  V  (t)  V  (l  + 

0 


r)  dt 


(A-6) 


Substituting  Eq.  (A-6)  into  Eq.  (A- 5)  gives 


(A-7) 


Equation  (A-7)  can  be  integrated  in  parts,  using  Leibnitz'  rule,  to  give 


(A-8) 


This  relation  is  given  by  Hinze  (Ref.  37)  who  credits  Taylor  (Ref.  66), 
and  Kampe  de  Feriet  (Ref.  67).  The  relation  is  also  given  by  Tennekes 
and  Lumley  (Ref.  68). 
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APPENDIX  B 
MIELV  PROGRAM 

A  program  (MIELV)  was  written  to  determine  the  characteristics  of 
laser  velocimeter  (LV)  signals  and  photomultiplier  tube  (PMT)  signal-to- 
noise  (S/N)  ratio.  A  general  Mie  scattering  routine  (SUBROUTINE  CALC) 
provided  the  Mie  scattering  functions.  In  particular,  MIELV  computes 
(1)  the  Mie  scattering  characteristics  of  spherical  particles  traversing 
a  plane  electromagnetic  wave  in  an  isotropic  medium,  (2)  the  scattered 
energy  incident  on  a  collector  lens  located  in  a  spherical  coordinate  system 
with  the  origin  at  the  particle  center,  (3)  the  amplification  and  atten¬ 
uation  of  the  signal  by  the  PMT  amplifiers  and  laser  line  filters  typical 
of  LV  systems  used  at  AEDC,  and  (4)  the  S/N  ratio  and  visibility  of  the 
resulting  LV  signals. 

It  is  assumed  that  the  Mie  scattering  characteristics  of  a  dual¬ 
scatter  LV  can  be  represented  by  the  scattering  from  an  equivalent  single 
beam  of  uniform  intensity.  The  Probe  Volume  Intensity  Program  (PVIP,  Ap¬ 
pendix  C)  was  developed  to  estimate  the  equivalent  intensities.  No  re¬ 
strictions  are  placed  on  particle  or  medium  refractive  indices  or  collector 
optics  viewing  angle.  It  is  further  assumed  that  the  amplification  of  the 
LV  signal  by  the  LV  electronics  is  a  linear  process.  Both  the  single-bean 
model  and  the  S/N  ratio  calculations  assume  that  the  half-angle  between 
the  intersecting  beams  foi'ming  the  LV  probe  volume  is  less  than  3  deg. 

In  addition,  only  flat  lenses  are  considered  and  the  product  of  re¬ 
fractive  index  and  particle  size  parameter,  8,  is  restricted  to  values  less 
than  1 ,000 . 


MIELV  is  based  on  the  theoretical  considerations  described  in 
detail  in  Section  3.0,  with  the  exception  of  the  integration  of  the  scat¬ 
tering  intensity  over  the  system  collector  lens.  This  integration  is 
performed  over  a  disk  of  radius  R  (representing  the  lens)  by  applying 
the  trapezoidal  rule  to  a  grid  network  of  11  x  11  points.  The  grid  is 
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defined  in  a  Cartesian  coordinate  system  (x,y,z)  with  its  origin  at  the 
center  of  the  lens,  and  the  grid  spacing  is  set  to  AX  =  R/5  to  define 
grid  points 


x  =—  R  +  (i  -  1)  AX 

ij 

=  -»+  (j  -  1)  AX 


(B~1 ) 


If  a  grid  point  falls  outside  the  disk,  new  coordinates  for  that  point 
are  defined  by 


‘j 


:  2  £  i  £  10  and  j  >  5 
:  2  £  i  £  10  and  j  £  5 


(B-2) 


:  i  =  1  or  11  and  j  >  5 

:  i  =  1  or  11  and  j  £  5 

Each  grid  point  on  the  lens  is  converted  from  the  Cartesian  coordinate 
system  (x,y,z),  with  origin  at  the  center  of  the  lens,  to  a  spherical 
coordinate  system  (r,0,<j>)  with  origin  at  the  particle  center,  as  pic¬ 
tured  in  Fig.  B-1.  The  lens  is  assumed  to  be  tangent  to  the  sphere  that 
passes  through  The  coordinates  of  each  grid  point  are  determined 

with  the  relations 


\*  =  L  sin  0  -  y  COS  6 


y*  -  X 

7*  =  L  cos  0  -  y  sin  8 


r 


+ 


*  2 


+  2 
y 


6* 


(B-3) 


<y>*  =  cf>  +  tan  ' 
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Z 


Figure  B-1.  Scattering  coordinate  system. 

With  the  point  (r*, ,  evaluation  of  functions  and  f^  is 
performed,  where  and  f9  are  defined  by 

f.  =  i  ■sin~  <£>*  _L 

f2  -  i n  fcin“  <f> ’  (B-4) 

After  evaluating  f^  and  at  a11  121  grid  points,  the  program  next 
does  the  trapezoidal  integration  by  summing  up  the  volumes  of  elements 
as  shown  in  Fig.  B-2. 


The  element  volumes  are  given  by  the  relation 

v  =  dx  dy  (fu  -  f,_,  f  j  +  fs_  i  .  ^  f,.  j-P 

where 


2d\  =  (x. 

h 

~  Vl'  j 

x i’  rl  Vi-1T  j-l') 

2dv  =  (y 
'  -  ij 

“  V'.  J-l 

1 

1 

[ 

1 

4 

(B-5) 


(B-6) 
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Figure  B-2.  Trapezoidal  integration  volumes. 

All  the  volumes  are  summed  to  approximate  the  integrals  over  the 
solid  collection  angle: 


fL 

sin^  4>*  dw 

<B—  7) 

sin"  cf>*  dw 

(B— 8) 

where  the  i  and  i  |(  intensity  functions  are  evaluated  at  the  center  of 
the  lens.  These  integrals  are  used  to  evaluate  the  total  energy  scat¬ 
tered  from  the  particle  onto  the  collector  lens,  as 


where 


Ec~  l\  +  [2 

(B— 9) 

'■  ■  (  \)A 

it  d2  V  2n  104  /J 

P 

(B— 1 0 ) 

'2'  *  d” 

(B-1 1 ) 

The  total  energy  reaching  the  PMT  photocathode  is  then  evaluated  using 


F  =  E. 


(B-1 2) 
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and  the  PMT  anode  current  is  determined  from 


i.  -  HC  '  F 


(B-13) 


However,  if  i  is  greater  than  the  saturation  limit  of  the  PMT,  i^  is 
set  to  the  saturation  limit.  The  total  signal  amplitude  is  given  as 

ed  ~  'a  Re  CB-14) 

Next,  the  theoretical  signal-to-noise  ratio  and  particle  velocity 
product  is  determined  from 


(S/N)  ■  V  =  Fj  F2  ]\  (rdp  x  10-4/K5)  (  - 


(B-1 5) 


where 


with 


Fj  =  V  G+  t/[(6.6  x  1034)  (3  x  1014)] 


f2  =  ec  Qc  A  Kj 


a  =  -si£ 

0  L2 


(B-1 6) 
(B— 17) 

(B-1 8) 


G*  =  G*  (G  -  1) 
C"  1C"'1'1  -  l) 


(B-1 9) 


and  where  =  i^  evaluated  at  the  lens  center  for  6=0,  whereas  k n 
i|l  evaluated  at  the  lens  center  for  S  =  0.  Finally,  the  theoretical 
particle  visibility  is  calculated  from 

J  j  (n-  x  HT4/ K*1 


FS  = 


frdpX  10 


(B-20) 


The  preceding  relations  form  the  basis  for  MIELV  which  was  pro¬ 
grammed  in  Fortran  IV  for  an  IBM  370/165  computer.  Detailed  documen¬ 
tation  on  the  program  is  available  upon  request  from  the  Commander, 
AEDC,  ATTN:  ACI. 
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APPENDIX  C 

PROBE  VOLUME  INTENSITY  PROGRAM  (PVIPJ 

This  program  obtains  the  local  illuminating  intensity  in  the  probe 
volume  of  a  dual-scatter  laser  velocimeter,  as  well  as  various  average 
intensity  values  for  the  fringe  model  of  Ref.  64.  According  to  the 
theory,  the  probe  volume  must  be  formed  by  the  intersection  of  two 
equal-intensity  laser  beams  each  having  a  Gaussian  light  distribution 
and  similar  polarization  and  intersecting  at  a  small  angle.  In  partic¬ 
ular,  the  probe  volume  is  assumed  to  be  an  ellipsoid  as  in  Fig.  C-1. 


Figure  C-1.  Ellipsoid  probe  volume. 

The  volume  of  this  ellipsoid  is  known  analytically  (V  =  ^  Trab2) ,  and 
the  accuracy  of  the  numerical  integration  can  be  checked  by  comparing 


to  the  analytical  result. 

The  average  intensity  in  the  probe  volume,  as  given  in  Section 
3.1.2,  Eq.  (77),  is 
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? '  Iff1'  ■  V  ■  t)  dx'  dv'  dz/< 


(C-2) 


where  _  2  2  2 

-2[X*  +  X*  -r  Z“  ] 

I1  =  I,  [0  =  [co^h  (4Y*Z*)  +  COS  (N>Y*)J  e  (C-3) 

''■  il li  X^I  =  X  |^/bt  Y  ,  =  Y  ^/h,  ami  Z  t  =  Z^a 

Also,  the  most  probable  intensity  is  given  in  Section  3.1.2,  Eq.  (79), 


where 


with 


<i>  jjifiohiw, 

,z 


AZ, 


i*  ^  _  ULl 


f  f  dxl  dYl 

*/y  */y 

1  1  Al 


(C-4) 


(C-5) 


(C-6) 


The  program  calculates  illuminating  intensity  at  various  Z* 
locations  as  Y*  varies  along  the  probe  volume  axis  from  Eq.  (C-3) . 
Then  the  average  intensity  is  calculated  as  it  varies  along  Z*,  from 
Eq.  (C-3). 


All  integration  is  accomplished  by  means  of  16-point  Gaussian 
quadrature.  The  triple  integration  is  performed  by  a  16-point  Gaussian 
quadrature  embedded  within  a  16-point  Gaussian  quadrature  which  itself 
is  embedded  within  a  16— point  Gaussian  quadrature.  Because  of  the  oscil¬ 
latory  nature  of  I,  the  integration  over  Y  is  divided  into  M  equal 
intervals,  each  of  which  is  performed  by  a  16-point  Gaussian  quadrature 
where 
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M  «  EN  +  l  (C-7) 

7 

This  formula  was  derived  after  examining  results  over  a  range  of  NN 
values  from  1  to  50. 

The  preceding  relations  are  incorporated  in  PVIP  which  was  pro¬ 
grammed  in  Fortran  IV  for  an  IBM  370/165  computer.  Detailed  computer 
program  documentation  is  available,  upon  request,  from  the  Commander, 
AEDC,  ATTN:  ACD. 


136 


AEDC-TH-78-65 


NOMENCLATURE 

A  Cross-sectional  area  of  probe  volume 

A,B,C  Constants,  Eq.  (37) 

A(X),  B(X),  C(X)  Constants,  Eq.  (39) 

A^,  Scattering  amplitude  coefficients,  Eqs.  (62) 

and  (63) 


a 


b 

b  (X) 


c 

P 


D 

d 

P 

E 

E 


e 


e 


d 


e 

s 


Semiminor  axis  of  probe  volume,  Eq.  (75) 

Waist  radius  of  laser  beam 

Mixing  layer  half  width,  Eq.  (38) 

Fluid  specific  heat  at  constant  pressure 

Jet  exit  diameter 

Particle  diameter 

Electric  field  vector 

Rate  energy  is  received  by  collector,  Eq.  (68) 

Exponential  (2.71828...),  or  hot-wire  output, 
Eq.  (37) 

Detected  voltage,  Eq.  (82) 

Supply  voltage  to  PMT 
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F<e,<t) 


f 

f(  ) 

<f> 


f 


D 


f 

o 


G. 

i 


G* 


g 


c 


h 

I 

I 

o 


Functions,  Appendix  B 

Angular  scattering  amplitude  function 

Measured  velocity  moment,  Eq.  (1)  or  lens  parameter 

Probability  density  function,  Eq.  (78) 

Arithmetic  average  velocity  moment,  Eq.  (3) 
Frequency  of  particle  fringe  crossings 
Bragg  system  reference  frequency 
Functions,  Eq .  (B-4) 

Stage  amplification 
Gain  relation,  Eq.  (B-19) 

Gravitational  acceleration 

Magnetic  field  vector 

PMT  transfer  function,  Eq.  (85) 

Constant,  Eq .  (37) 

Illuminating  light  intensity 

Peak  single  beam  light  intensity,  Eq.  (73) 

Average  probe  volume  light  intensity,  Eq.  (77) 
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I* 


I* 


<1> 

<1  > 
max 

ij_(e),  i|((0) 


J 

j 

c 

J 

n 

K 


L 


T 


Nondimensional  illuminating  light  intensity,  Eq.  (76) 

Average  intensity  along  the  constant  Z*  planes, 

Eq.  (80) 

Most  probable  probe  volume  intensity,  Eq.  (79) 

Most  probable  maximum  probe  volume  intensity,  Eq.  (81) 

Scattering  intensity  functions,  Eqs.  (66) 
and  (67) 

LV  signal  current,  Eq.  (83) 

Unit  vectors 

Scattered  light  intensity,  Eq.  (55) 

Energy  unit  conversion  from  Btu  to  ft-lbf,  Eq.  (36) 

Bessel  function  of  the  first  kind 

Kurtosis,  Eq.  (35) 

Fringe  spacing 

Processing  delay  time 

Distance  from  probe  volume  to  collector  lens 
Integral  time  scale,  Eq.  (13) 
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M 


m 


m 


1 


m 


2 


N 

NN 


n 


P 


P<  ) 


P  (  ) 
n 


R 

R(t) 

R 

e 


Number  of  samples  in  a  data  set  or  stages  of 
amplification  or  integer,  Eq.  (C— 10) 

Refractive  index  ratio 

Particle  refractive  index 

Medium  refractive  index 

Particle  number  density  or  stages  of  attenuation 
Integer,  Appendix  C 

Number  of  fringes  in  the  Z*  =  Oi probe  volume  plane 
Summation  index 

Pressure,  Eq.  (36),  or  point  location.  Fig.  25 
Probability 

Laser  power 
Legendre  polynomials 
Riccati-Bessel  function,  Eq.  (64) 

Collector  lens  radius 

Velocity  autocorrelation  coefficient,  Eq.  (11) 

LV  circuit  resistance 
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fAt 


RNVn 


^VAt 


^VAt* 


Correlation  coefficient,  Eq.  (7) 

Correlation  coefficient,  Eq.  (53) 

Correlation  coefficient,  Eq.  (43) 

Correlation  coefficient  for  processing  delays 


r 


roOO 


r 

a 


S 

S 


s 

c 


s 


T 


Radial  dimension 

Radius  of  jet  core 

Particle  radius 

Core  radius  for  turbulence  intensity  curve  fit, 
Eq.  (39) 

Skew  coefficient,  Eq.  (23) 

Poynting  flux  vector 

Collector  lens  surface  area 

Number  of  realizations  in  an  interval,  Eq.  (49) 
Number  of  samples 

Complex  amplitude  functions,  Eqs.  (58)  and  (59) 
Variable  of  integration,  Appendix  A 
Time  interval,  Eq.  (1) 
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T. 

1 


u,v,w 


V 


A 

V. 

1 


p 

|v| 

<v> 


A 

V. 

1 


V  V 

Vvi 


X 


WZ1 


Stage  attenuation 
Sampling  time 
Temperature 

Time 

Initial  time 

Velocity  in  the  main  flow  direction 
Jet  exit  velocity 
Fluctuating  velocity  components 
Velocity 

Temporal  mean  velocity  over  interval  At^ 
Volume  of  probe  volume  region 
Velocity  magnitude,  Eq.  (25) 

Arithmetic  average  velocity 

Mean  velocity  over  the  interval  At^ 

Sinusoidal  velocity  parameters,  Fig.  20b 

Axial  location  from  nozzle  exit 

Probe  volume  coordinate  system.  Fig.  33 
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X*  Y*  Z* 

Vi  *  i 

x,y,z 

x* 

a,B,Y 

a 

a . 

1 

y 

AP 

At 

At* 

<At> 

Au. 

l 

A<V> 

AX 

6N 

n 


Nondimens ional  probe  volume  coordinates 
Cartesian  coordinate  system,  Fig.  25 

Function 

Error  coefficients,  Eqs.  (30-32) 

Particle  size  parameter 
Velocity  parameter,  Eq.  (27) 

Particle  characterizing  parameter 
Ratio  of  specific  heats,  Eq.  (36) 

Differential  pressure,  Eq.  (36) 

Time  interval  between  realizations 

Delay  time  after  acquiring  a  data  point,  Eq.  (40) 

Arithmetic  average  time  interval 

Velocity  precision,  Eq.  (29) 

Confidence  interval  for  the  mean,  Eq.  (33) 

Finite  step  size,  Appendix  B 

Number  of  particles  entering  probe  volume,  Eq.  (8) 
Fraction  of  laser  power  reaching  the  probe  volume 
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e 

x 

\ 

n 


K 


TT 


¥ 

n 


P 

E 


CT- 

V 


0 


<v> 


A  o 


<v> 


T 


X 

n 


$ 

n 


$ 


Convergence  angle  of  laser  beams.  Fig.  1,  or  scattering 
angle.  Fig.  25 

Wavelength 

Sampling  rate 

Series  multiplication 

Irrational  number  (3.141592...) 

Angular  dependence  function,  Eq.  (60) 

Fluid  density 

Summation 

Time-integrated  turbulence  intensity,  Eq.  (21) 

Arithmetic  average  turbulence  intensity 

Confidence  interval  for  the  intensity,  Eq.  (34) 

Time  difference,  Eq.  (10),  or  Doppler  burst  period, 

Eqs.  (30)- (32) 

Angular  dependence  function,  Eq.  (61) 

Riccati-Bessel  function,  Eq.  (65) 

Phase  angle.  Fig.  20b,  or  scattering  coordinate, 

Fig.  25 
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Collector  solid  angle 
Frequency,  Fig,  20b 


Ambient 


Centerline 


Cell  gage  conditions 

Indices 


Maximum 

Summation  indices,  Eq.  (18) 
Perpendicular  component 

Parallel  component 


SUPERSCRIPTS 

Approximate-time-integrated  average,  Eq.  (2) 
A  True  time- integrated  average,  Eq.  (1) 

■*  Fluctuating  component 
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